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ABSTRACT 



We investigate the dynamical evolution of the Galactic Globular Cluster 

^ ■ System in considerably greater detail than has been done hitherto, finding that 

^ ■ destruction rates are significantly larger than given by previous estimates. The 

O ' general scheme (but not the detailed implementation) follows [Aguilar, Hut, fc 

g : I Ostrike^ (1988; AHO). 

^ ■ For the evolution of individual clusters we use a Fokker-Planck code 

■ including the most important physical processes governing the evolution: 

p two-body relaxation, tidal truncation of clusters, compressive gravitational 

Q ■ shocks while clusters pass through the Galactic disk, and tidal shocks due to 

^ . passage close to the bulge. Gravitational shocks are treated comprehensively, 

ci ■ using a recent result by |Kundic fc Ostriker| (1995) that the {AE"^) shock-induced 



relaxation term, driving an additional dispersion of energies, is generally more 
^ . important than the usual energy shift term (AE). Various functional forms 

. of the correction factor are adopted to allow for the adiabatic conservation of 

stellar actions in a presence of transient gravitational perturbation. 

We use a recent compilation of the globular cluster positional and structural 
parameters, and a collection of radial velocity measurements. Two transverse to 
the line-of-sight velocity components were assigned randomly according to the 
two kinematic models for the cluster system (following the method of AHO): one 
with an isotropic peculiar velocity distribution, corresponding to the present day 
cluster population, and the other with the radially-preferred peculiar velocities, 
similar to those of the stellar halo. We use the Ostriker & Caldwell (1983) and 
the Bahcall, Schmidt, & Soneira (1983) models for our Galaxy. 

For each cluster in our sample we calculated its orbits over a Hubble time, 
starting from the present observed positions and assumed velocities. Medians of 
the resulting set of peri- and apogalactic distances and velocities are used then 
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as an input for the Fokker-Planck code. Evolution of the cluster is followed 
up to its total dissolution due to a coherent action of all of the destruction 
mechanisms. The rate of destruction is then obtained as a median over all the 
cluster sample, in accord with AHO. 

We find that the total destruction rate is much larger than that given by 
AHO with more than half of the present clusters (52% — 58% for the OC 
model, and 75% — 86% for the BSS model) destroyed in the next Hubble time. 
Alternatively put, the typical time to destruction is comparable to the typical 
age, a results that would follow from (but is not required by) an initially power 
law distribution of destruction times. We discuss some implications for a past 
history of the Globular Cluster System, and the initial distribution of the 
destruction times raising the possibihty that the current population is but a 
very small fraction of the initial population with the remnants of the destroyed 
clusters constituting presently a large fraction of the spheroid (bulge+halo) 
stellar population. 

Subject headings: celestial mechanics, stellar dynamics - globular clusters: 
general - Galaxy: kinematics and dynamics 

1. Introduction 

Globular clusters are thought to be the oldest stellar systems in our Galaxy, and a 
history of attempts by theoreticians and observers to understand the keys of their evolution 
is as old as the discipline of stellar dynamics, with a comprehensive review provided by 
Spitzer] (1987). Yet, these relatively simple stellar systems are not understood well enough 
to predict their future with desirable accuracy, in part due to lack of accurate observational 
values for the current dynamical state and in part due to a residual uncertainty concerning 
the complete catalog of relevant physical processes operating on these systems. More than 
that, we have almost no clues about their past, and in particular, whether what we see 
now is representative of the initial Globular Cluster System, or just a small leftover after 
a great destruction battle that occurred earlier in the history of the Galaxy. In this paper 
we consider the evolution of the globular clusters and their ultimate disappearance, and 
propose a simple model for their initial distribution. 



Pre- and post-core collapse evolution of an isolated cluster is relatively well understood 
( [Spitzer] 1987; [Goodman] 1993). Significant progress in understanding of the evolution was 
achieved using Monte Carlo and Fokker-Planck calculations. But the galactic environment 
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makes the clusters subject to external perturbations - tidal truncation and the gravitational 
shocks due to passages close to the bulge and through the disk. The shock processes, 
although known to be important, have never been carefully included in the evolution of the 
system. We investigate the Fokker-Planck models including the shocks elsewhere (Gnedin 



Lee, fc Ostriker 1996| ; hereafter GLO) and show that the dispersion of energy of the stars. 



induced by the shocks ( [Kundic fc Ostriker] 1995; hereafter KO), is generally even more 



important for the evolution than the first-order energy shift. Another very important aspect 
of modelling globular clusters is the initial mass spectrum. Multi-mass clusters undergo 
core collapse much faster than in the single-mass case, and their destruction is much more 
efficient (see |Lee fc Goodman] 1995, and the references in GLO). We restrict ourselves 



however to the single mass models in order to maintain clear physical understanding; but 
due to omission of this aspect of the problem, our results provide a lower bound to the rate 
of destruction of the globular clusters. 

As long as we have (at least, approximate) understanding of the evolution of a 
single cluster we can turn to the study of the system of globular cluster in our Galaxy, 
and in external galaxies. ]Ghernofl:', Kochanek, fc Shapiro (1986)] used a semi-analytical 



Monte-Carlo technique to estimate the importance of the different mechanisms acting upon 
the cluster. They considered two-body relaxation, tidal stripping of stars, and the first order 
tidal shocking effect - due to the crossing disk and interactions with giant molecular clouds 
(GMC). For each of those processes they calculated the cluster mass and energy changes 
associated with them to predict the evolution. They followed the cluster evolution only up 
to core collapse, and assumed a single-mass King model (]King 1966|) for the internal cluster 
structure. Tidal heating due to the GMC was found to be negligible compared to the disk 
shocks. Note however, that the Galactic model assumed in that work ( Pahcall, Schmidt, fc 



Soneira 1983] ; BSS) is strongly favorable to the disk shock for the clusters with small orbital 



radius since the surface density of the disk increases exponentially as the galactocentric 
radius decreases (see Section p. 2] ). Chernoff et al. (1986) has concluded that many of the 
clusters located within inner 3 kpc from the Galactic center have undergone core collapse, 
and many of them may already have been destroyed. A number of authors have pointed 
out that the bulge and stellar spheroid themselves could be composed of remnants of the 
destroyed globular clusters, if prior destruction of globular clusters occurred at high enough 
rate. 



Another mechanism for the mass loss is the stellar evolution, phernoff' fc Shapiro 



(1987) used a similar method to Chernoff et al. (1986) and included a power-law initial 
mass function for the cluster stars. Mass loss due to stellar evolution is important for the 
early evolution of the cluster, but then fades away because the mass loss is large only for 
massive stars whose lifetime is short. Comprehensive study including the effects of stellar 
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evolution has been done by Uhernotf fc Weinberg (1990)| . They used a Fokker-Planck code 



with an extensive spectrum of stellar masses (20 species), which includes stellar evolution 
and relaxation processes. Multi-mass models evolve much faster and the evaporation rate is 
larger. Also, mass segregation (|Spitzer| 1987) speeds up the collapse. Mass loss during first 



5 X 10^ yr is sufficiently strong to disrupt weakly concentrated clusters (c < 0.6). Combined 
with the relaxation, it destroys many low mass and low concentration clusters within a 
Hubble time. 

The present characteristics and the evolutionary state of the observed Galactic globular 
clusters were also investigated by [Aguilar, Hut, fc Ostriker] (1988; hereafter AHO). We will 
draw heavily on that paper and compare our results with AHO. AHO used a sample of 
the 83 Galactic globular clusters with the known structural parameters and line-of-sight 
velocities. They considered virtually all the important physical mechanisms (except 
the mass spectrum and mass loss due to stellar evolution) to calculate the present day 
destruction rates for the clusters in the sample. The rates were defined as the inverse time 
it takes for a given mechanism to dissolve a cluster, in units of a Hubble time, which is 
nominally defined as lO^'' yr. They estimated the rates for the evaporation through tidal 
boundary, disk and bulge tidal shocks, and dynamical friction. The last process has not 
been widely investigated in its application for the clusters' evolution. Its effect reduces to 
the gradual spirahng of the cluster toward the Galactic center, as the cluster loses orbital 
energy due to continuous interactions with field stars and dark halo. AHO found that 
this mechanism is a relatively unimportant one, except for unusually massive clusters. 
They calculated a number of orbits associated with each cluster in the sample and took a 
median over the whole resulting distribution for a corresponding destruction rate. They 
investigated two Galactic models (OC and BSS; Section |2?^ ), and two kinematic models 
for the globular cluster system: isotropic, which corresponds to the current distribution, 
and predominantly radial, which resembles that of the halo stars and might be closer to an 
initial cluster distribution. The assumption concerning the velocity distribution is necessary 
as only one velocity component is observed. (Note however that the program to obtain the 
true space velocities of globular clusters is underway; see [Gudworth fc Hanson] 1993). We 
will discuss the kinematic models in Section |2.3| . 

The central bulge was found to be very efficient in destroying clusters on a highly 
elongated orbits. This leads to an isotropization of the orbits, and even preferential survival 
of tangentially biased orbits in the Hubble time. Thus, an initially radial distribution better 
fits the present population, than the initially isotropic one. At the current time, evaporation 
is the most important destruction process. The difference between the two Galactic models 
was found to be relatively small, except for the disk shocking that obviously depends on 
a body of the surface mass profile. AHO introduced a weighting factor that accounts for 
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the fact that some orbits are less probable because of the various destruction mechanisms 
acting upon the cluster. Such a weighting reduces the computed rates, but converges to a 
similar value of 0.05 for all Galactic models and kinematic profiles considered. This implies 
that only 4 or 5 clusters are destroyed over a Hubble time for present conditions. But AHO 
treated clusters in a simplified fashion as King models and they (like other investigators) 
did not allow for the tidal shock relaxation phenomena. It is these defects that we remedy 
in the present paper. 



Finally, a recent paper by [Hut fc Djorgovski (1992)] used a statistical approach and 



approximate analytical estimates for the relaxation times for a sample of 140 clusters. They 
concluded that the current evaporation rate is 5 ± 3 Gyr~^ (about ten times the rate found 
by AHO), which means that a significant fraction of the present day clusters would be 
destroyed in a next Hubble time. 

In this paper we make a significant improvement over AHO by applying detailed 
Fokker-Planck calculations to the real globular cluster sample. We describe the sample in 



Section p.l| and the Galactic models in Section p.2| . Then we discuss the two kinematic 
models and compare the resulting properties of the orbits for our sample. In section |2.5| we 
describe the destruction mechanisms involved in the globular cluster evolution. We conclude 
section § with the history of our code and the formulation of the numerical strategy. Section 
H presents our results for all runs. Finally, we speculate on the possible past history and 
future fate of the Galactic globular clusters in Section ^. Section || sums up our conclusions. 



2. Method 
2.1. Observational Data 

We have used a recent compilation of globular cluster coordinates by PJorgovski fc 



Meylan (1993)| , and distances and absolute magnitudes from PjorgovsTci] (1993). Cluster 



concentrations, core radii, and half-mass radii were taken from [Irager, Djorgovski, fc King 



1993 )[ . The core-collapsed clusters were assigned a limiting concentration of c = 2.5. We 
used a constant mass-to-light ratio to obtain the cluster masses, {M/L)v = 3 in solar units 
( Pjorgovskil 1993). 

We selected 119 clusters out of 143 with available photometric data. All clusters in our 
sample have measured radial velocities which we collected from several sources. Pudwortli 



|fc Hanson] (1993) have derived actual space velocities for 14 clusters, but we chose not to 
use those data to maintain homogeneity of the sample. We assign two unknown velocity 



components to each cluster using our statistical method in Section pT3. 
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The observed parameters of our sample are given in Table 1. The first two columns 
are the sequential number and cluster's name. The next two columns are the Galactic 
coordinates. D is the distance in kiloparsecs of the cluster from the Sun, and R is its 
galactocentric radius. The seventh column is concentration c = log{Rt/Rc), followed by core 
radius Rc in parsecs. Tidal radius Rt was calculated from the previous two quantities. The 
next column is the cluster mass converted from its luminosity. The last two columns are 
the line-of-sight velocities and the reference numbers to the various sources. The references 
are listed in the comments to the Table 1. 



2.2. Galactic Model 

To evaluate the gravitational shocks on the clusters we have used two models for our 
Galaxy: the Ostriker-Caldwell (1983) model D-150 (OC), and the Bahcall, Schmidt, & 
Soneira model (1983; BSS). AHO present a detailed comparison between the two models. In 
the OC model the disk is represented as the difference of two exponentials which vanishes at 
the center, whereas BSS's disk density rises monotonically as R decreases. The BSS model 
has also a compact nuclear component within 1 central kiloparsec. Thus both disk and 
bulge shocks (Section p75| ) are expected to be more prominent for BSS model. As suggested 
by AHO, we consider the difference of the results or the two models as a measure of the 
uncertainty associated with the distribution of mass within the Galaxy. 

We use the same routine to compute the Galactic model as was done by AHO (the 
routine was kindly provided to us by L. Aguilar). 



2.3. Velocity Distribution for the Globular Cluster System 

Kinematic models for the Galactic Globular Cluster System (GCS) have been sought 
widely in the past (e.g., |brenk fc WEit^ 1980; [l'homas| 1989). |brenk fc WEit^ (1980) 



studied a sample of 66 clusters and found no evidence for the radial expansion of the 
GCS as a whole. Their best fit to the observed kinematic distribution is that with the 
isotropic velocity dispersion rising with the galactocentric distance as R^'^. An isotropic and 
isothermal (with constant velocity dispersion) model is still consistent with their data at the 
90% confidence level. Later work by |Thomas| (1989) included 115 clusters and confirmed 
the absence of the expansion. He found that within the inner 7 kpc to the Galactic center 
the velocity distribution is isotropic, but for the outer radii a velocity ellipsoid with slightly 
increasing line-of-sight velocity dispersion may be preferred. Rotation velocity estimates 
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and line-of-sight velocity dispersions as well as their errors for the both sources are given in 
Table |. 

We repeated the analysis of the velocity distribution using our sample. We adopted the 
solar motion relative to the LSR of (—9, 12, 7) km s~^ ([Mihalas fc Binney 1981| ), and the 
circular velocity of the LSR of 220 km s^^. Line of sight projections of the two velocities 
were subtracted from the observed radial velocities to obtain ones relative to an inertial 
frame. Then the rotation velocity and its uncertainty were estimated using equations (15) 
and (16) of [Thoma4 After subtraction of the rotation velocity along the line of sight we 



end up with the peculiar velocities with (presumably) zero mean. Both our v^ot and aios are 
in agreement with the above results within errors. Our adopted values are also summarized 
in Table ^ 

Following AHO we investigate two initial kinematic models for the globular cluster 
system. The first one has an isotropic peculiar velocity distribution with constant velocity 
dispersion. This is the simplest model still consistent with the observations. We used the 
one-dimensional dispersion of 118 km s~^, resulting from our sample. 

The second model is anisotropic with velocity ellipsoid axis ratio at solar position 
Tq = 8.5 kpc similar to the spheroid stellar population II: 

where ar and at are the one-dimensional dispersions radial and transverse to the 
Galactic center, respectively. We have adopted the AHO value of the anisotropy radius 
Ta = O.8^/^r0 7.6 kpc. Cluster orbits for this distribution are nearly isotropic within r^, 
and become more and more radial with increasing distance from the center. 

The amplitude of the radial velocity dispersion can be obtained from Jeans equation 
for the constant circular velocity potential. Integrating the Jeans equation in spherical 
coordinates gives (e.g., [Ogorodnikov 19651) : 

For the velocity ellipsoid given by equation ([l|) and the power-law density profile p oc r~", 
we obtain 

^2 /„,2 „,2 \ a~2 ~'~ a /o\ 

= (^circ - ^rot) „2 i 2 " 



a 



We accept f circ = 220 km s~^ as the standard value for the circular velocity. [Xhomas) found 
a = 3.5 for his sample of globular clusters. In this paper we assume a = 3, corresponding 
to the old spheroid stellar population. 
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The distribution of peculiar velocities at a given galactocentric position r of a cluster 
is then (AHO eq. [4]) 



fivr,vn,vt2) = A exp 



(4) 



2cr? 

where Vr is the radial galactocentric velocity, and fti, Vt2 are the two transverse components. 

The distribution of orbital eccentricities is shown in Figures |I| and |^ for the OC and 
BSS models, respectively. Upper panels in both figures show histograms for the anisotropic 
velocity ellipsoids, and the lower ones for the isotropic distribution. Note how the number 
of large eccentricities changes from the former to the latter. 



2.4. Orbit Integration 

In order to model the gravitational shock experienced by a cluster flying by the galactic 
center (bulge shock) we need an estimate of a perigalacticon distance Rperi and orbit shape 
(eccentricity e) for each cluster from our sample. We use the orbit integration routine 
described in AHO, which employs a 4th order Runge-Kutta method with variable time 
step. The Galactic model extends up to a distance of 250 kpc, and all orbits beyond that 
point are discarded. 

Only 4 out of the needed 6 phase space parameters required to start the calculation 
of orbits are known from the observations (Table 1): 3 spatial coordinates and a velocity 
component along the line of sight. To complement the two "missing" tangential velocities 
we used a statistical approach similar to that of AHO. The two velocities were randomly 
drawn and then selected using the rejection method ( [Press et al. 19921) in order for the 



total three-dimensional peculiar velocity be consistent with the assumed kinematic model 
(section |2.3| ). The systemic rotation velocity was then added to the peculiar velocity, and 
the resulting velocity with respect to the Galactic center was used for the orbit integration. 

We have followed trajectories that each cluster makes in 10 billion years. Perigalactic 
and apogalactic distances were calculated as medians from the set of all orbits. Similarly 
we determine the velocity of the cluster at the perigalacticon (which we use to calculate 
the amplitude of the bulge shock, section |2.5.3| ), and the vertical velocity component at the 



point where the cluster crosses the disk (to estimate the strength of the disk shock, section 



2.5.2). Thus only one "mean" orbit is used to evaluate the gravitational shocks. 



In contrast to AHO, we could not follow every orbit of a cluster because of the 
nature of the Fokker-Planck calculations. The code has the time step controlled by the 
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relaxation processes, and it seems hard to reconcile with the orbit integration procedure. 
We plan however to return to this subject in next paper and try to perform simultaneous 
Fokker-Planck and orbital calculations which would allow us to model the evolution of the 
clusters in a more natural way. 



2.5. Dynamical Processes 

2.5.1. Evaporation 

Two-body relaxation leads to the escape of stars approaching the unbound tail of the 
cluster velocity distribution ([Ambartsumian 1938| ; ppitzer 1940|) . Tidal truncation due to 



the Galactic potential accelerates this process. Much more dramatic effects result from 
the gravothermal instability, when the inner part of cluster contracts (core collapse), and 
the envelope expands. This, in turn, accelerates the rate of evaporation of stars from the 
cluster. A recent review of pre- and post-core-collapse evolution of a tidally truncated 
cluster is given by [Goodman] (1993). 

It is conventional to express the life-time of a cluster in terms of the half-mass 
relaxation time ([Spitzer fc Hart 1971|) 



where M is the total cluster mass, Rh is the half-mass radius, m* - average stellar mass, 
and ln(A) = ln(0.4A^) - Coulomb logarithm, N being the number of stars in the cluster. 
Henon 196T| introduced the "escape probability" 



c dM 

= -M -IT- 

Here dM/dt is the mass-loss due to two-body relaxation. For the self-similar solutions 
Henon found ~ 0.045. See GLO for detailed discussions of the evaporation time scale. 



2.5.2. Disk Shock 

When a cluster passes through the Galactic disk, it experiences a time-varying 
gravitational force pulling the cluster toward the equatorial plane. The characteristic 
time-scale for this force is the time it takes the cluster to cross two vertical scale-heights of 
the disk, H: 

Across — -r r ; (7) 
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and Vz is the cluster velocity component perpendicular to the Galactic plane. Because 
of the large orbital velocity of the cluster and the relatively small vertical extent of the 
disk, the crossing time is usually shorter than the orbital period of stars in the outer part 
of the cluster. Owing to the short-term nature of the effect, it was called a "compressive 
gravitational shock" QOstriker, Spitzer, fc Chevalier 19721) . On average, stars gain energy 



and the cluster binding energy is reduced. This accelerates the escape of stars from the 
cluster through evaporation. 

On the other hand, in the central region of the cluster the effect of the compressive 
shock is largely damped, since the stars move very fast and their orbits become adiabatically 
invariant. The impact of the shock is thus a strong function of the position of a star inside 
the cluster. Following |Spitzer|, we define an adiabatic parameter 



2ujH 



where uj is the angular velocity of stars inside cluster (assuming for simplicity circular 
motions), and the subscript d stands for disk; Xd represents the ratio of the shock duration 
to the orbital period of a star, so that for small values of uj the term "shock" is appropriate. 

To evaluate quantitatively the effect of the disk shock, several approaches have been 
proposed in the literature. The simplest one, the impulse approximation, assumes that the 
shock is so fast that the star does not change its position in the cluster significantly over 
the time t cross- However, for the reason of adiabatic conservation, this approximation highly 
overestimates the impact when Xd ^ 1- A more careful treatment is given by a harmonic 
potential approximation ([Spitzer 1958| ), where we assume all stars, initially at same radial 



distance r from the center of the cluster, move around the center with the same oscillation 
frequency uj. Then, referring to the equation (25) of KO, the average energy shift for every 
star is 

(Ai?)disk = ^||^A,(x,), (9) 

where Qm is the maximum gravitational acceleration experienced by stars due to the disk. 
Here Ad{xd) is the factor taking into account the adiabatic invariants. For the harmonic 
approximation (e.g. |Spitze4 eq. [5-28]) 



^,i(x)=exp -y . (10) 



We refer to equation (^) as a Spitzer correction. However recently Weinberg (1994) 



showed that small perturbations of stellar orbits can still grow in a system with more than 
one degree of freedom. If the system is represented as a combination of multidimensional 
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nonlinear oscillators, it is very likely that some of the perturbation frequencies will be 
commensurable with the oscillation frequencies of stars. Then those orbits receive a 
significant kick from the perturbation, and thus no longer conserve their actions. Averaging 
over whole cluster can give an appreciable change in velocity and energy. Such resonances 
can occur even for an arbitrary small resonant mode. Consequently, the adiabatic factor 
Ad{x) is not exponentially small for large x, but rather a power-law. The simplest form of 
the correction can be written as 

/ 2 \ -3/2 

A,ix)=^l + ^j . (11) 

For large x ^ 1, oc x~^.^ In the following, we call equation (|ll]) the Weinberg correction. 
Figure |^ shows the difference between the two forms of the adiabatic correction. A much 
more detailed discussion of the Weinberg correction can be found in [Gnedin et al.| (1996). 



Besides the shift in energy of stars (AE), a gravitational shock also induces a quadratic 
term {AE"^), which governs the dispersion over energy spectrum and thus pushes some 
loosely bound stars outside the tidal radius, speeding the disassociation of the cluster. 
"Shock-induced relaxation" was mentioned briefly in [Spitzer fc Chevalier] (1973; see also 



Spitzer 1987, p. 1160). Recently KO noted the importance of this effect. We refer a reader 
to that paper, in which it is shown that this tidal shock relaxation can be in many cases 
competitive with two-body relaxation in causing the evolution of a cluster. Like ordinary 
relaxation it causes a spread (diffusion) of initially similar orbits and ultimately will tend 
to induce core collapse. 



According to KO (their eq. [25]) 

_ 4/ 



(AE2),,3, = 14^^,2(0;), (12) 



where in the harmonic approximation 

Ad2{x) = - exp . (13) 

Note that the correction in this form does not match unity for a; <^ 1. Rather it goes a 
factor 1.8 over the impulse approximation. The harmonic potential obviously does not 
apply for the outer parts of the cluster, so some mismatch is expected. On the other hand. 



^This power law index was suggested by S. Tremaine. 

^We are indebted to L. Spitzer for pointing out this reference to us. 
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it enhances the shock in the cluster halo. We have chosen not to modify the Spitzer's 
formula, and compare the calculations with those resulting from the Weinberg correction. 
The latter has been used in the same form of equation (^) for both energy terms. We have 
performed calculations with the two forms of the adiabatic correction, and also a test case 
in the impulse regime (no correction was applied. Ad = 1). 

Disk shocks occur twice during the orbital period of the cluster, so we define the disk 
shock time scale at the half-mass radius (neglecting the adiabatic vcorrections) as 



. ^ PoA ( -E 




(14) 



where Porb is the orbital period, E ^ —0.2GM/Rh - cluster energy, and the AE and u are 
evaluated at Rh- Analogously, we find that 



Porb \ 3 



W,2 = ( ) ~ ^^disk- (15) 



Both these terms contribute to the destruction of the cluster (although by quite different 
processes: (AE) enhances evaporation and {{AE)"^) enhances core collapse), so that the 
total destruction rate associated with the disk shock may be written roughly as 

1 1 7 _i 
i/di,k = - — + = - tdisk- (16) 

^disk rdisk,2 ^ 



2.5.3. Bulge Shock 

Similar to the compressive disk shock, every globular cluster experiences a tidal shock 
during its passage close to the Galactic center. The massive compact component at the 
center of the Galaxy (bulge) induces a strong tidal force on the cluster near the perigalactic 
point of the cluster orbit. The difference between this effect and that from the smooth and 
steady tidal field of the Galaxy is primarily due to the time dependence of the bulge shock. 

The very close effect of the tidal shock induced by giant molecular clouds was 
considered by L. Spitzer as early as (1958). The disturbing object was represented as a 
point-mass, and cluster orbit near the point of the closest approach (perigalacticon in our 
case) is assumed to be a straight path. Employing the harmonic approximation we find the 
energy change for each star due to the bulge shock: 

(Ai5;)buige = ^ {y^^ ^w(x,) xiRp) KRp, Ra)- (i7) 
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Here is the bulge mass, Vp is the cluster velocity at the perigalacticon Rp, and Ah is the 
corresponding adiabatic correction. Two new corrections arise as follows. The distribution 
of the bulge mass extends up to many kiloparsecs, and for some clusters it is not a good 
approximation to consider the bulge as a point-mass. Rather, the tidal field exerted by 
such an extended mass profile will differ from the point-mass field (given by, e.g., ^pitzei) ) . 
Details of the calculation of the correction factor x{Rp) allowing for that effect will be given 
elsewhere ( pnedin fc Hernquist 199^ ). 

We give here only the final expression which we use in our calculations: 

x{Rp) = \ [(3Jo -h- hf + (2/o -h- 3 Jo + + /o'] , (18) 

where 

h{Rp) = / m6(i?/) — — , (19a) 



(C^ - 
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h{Rp) ^ MRpO ^2 ((2 _ 1)1/2 ' (19b) 
MRp) ^ / ruhiRpO .4 ..2 \m/2 ' (19c) 



MRp) ^ MRpO ^4 ((2 _ 1)1/2 > (19d) 

and mh{R) = Mf,{R)/Mh is the normalized bulge mass distribution at radius R, with 
mb{R) = dmb{R) / din R. 

Aguilar & White (1985; also AHO) proposed a correction of the form x = 72 -^o only. 
Obviously that correction factor is everywhere smaller than our's (eq. [0]), although the 
difference is not dramatic. We have plotted both these factors in Figure ^ For comparison, 
we plot also the function mb{R), which would be the naive correction which allowed only 
for the mass interior radius R. 

The second correction. A, is intended to take into account the time variation of the 
tidal force along an elliptic orbit of a cluster. Following AHO, we take the difference of the 
magnitudes of the tidal force at perigalacticon and apogalacticon, as the total amplitude of 
the tidal effect. Thus the correction factor A can be expressed as (AHO eq. [11]) 



\{Rp, Ra) 



Mb{Ra) [R, 



up 



(20) 



Mh{Rp) \Ra 

where Rp and Ra are the perigalactic and apogalactic distances of the cluster, respectively. 

Finally, Ahi{x) is the corresponding adiabatic correction for the bulge shock. In the 
harmonic approximation we get the Spitzer correction: Ah{x) = ^j^ \Lx[x) + Ly{x) + Lz{x)] 
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in his notation (|Spitzer 1958 , eqs. [36-38]). For the Weinberg correction we use again the 



same function (eq. |]rT|) as for the disk shock. The argument of the A function is now 
different, 

= (21) 

Figure ^ compares the two corrections. 

As in the case of the disk shock, bulge shock induces the second relaxation term 

8 f C M \ 

(Ai?')buigc = - (^^j u'r'A,^{x,)x{Rp)mp.Ra)- (22) 

Analogously to the disk case, Ai)2 = O/SAfei in the Spitzer regime, and ^4^2 = ^ti for the 
Weinberg's one. 

Similarly we define the bulge shock time scale. In this case however, the effect occurs 
only once per cluster orbital period. 

-E\ 3 (VpRl^'' 



ibulge,2 = -Porb (^^1 = ^ ^bulge- (24) 



Finally, total destruction rate associated with the bulge shock is 

1 1 7^-1 



t'bulge — + — ^ ^bulge- (25) 



■^bulge ^bulge,2 3 



2.6. Fokker-Planck code 

We calculate the dynamical evolution of the globular clusters from our sample using 
an orbit-averaged Fokker-Planck code descended from that of Cohn (1979, 1980). The code 
has been modified by [Lee fc Ostriker (1987)| and [Lee, I'ahlman, fc Richer (1991)| to include a 



tidal boundary and three-body binary heating. Although the code allows a multicomponent 
stellar mass function, we restricted ourselves in this paper to a single-component case 
[m^, = 0.7 ruQ). This requires less parametrization of the numerical models and, we hope, 
gives clearer understanding of the physical processes involved. In the future, it certainly 
will be of interest to generalize the present calculations, including a realistic mass function. 
Lee fc Goodman!, for example, considered evaporation of a multi-mass cluster in a steady 



tidal field and found that the mass loss doubles compare to the single-component case. 
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Stars beyond the tidal boundary are not lost instantaneously, but rather follow 
continuous distribution function f{E), as described in |Lee fc Ostriker (1987)| . This takes 



into account the fact that the tidal radius is not a strict "border" for the cluster, because 
the internal force just balances the Galactic tidal force at that point. Therefore the outer 
stars escape only when they go further away from the cluster. The tidal field of the 
Galaxy is assumed to be steady and spherically symmetric. The latter is a weakness of a 
one-dimensional code (for a discussion see [Lee fc Goodman] 1995). The former assumption 
is valid only for a circular cluster orbit. For the actual elliptical orbit the maximum tidal 
stress occurs close to the perigalactic point, and probably determines the tidal cutoff radius. 
A theoretical estimate for Rt is given by [Innanen, Harris, fc Webbink (1983)| who assumed 



a spherical mass distribution that increases linearly with galactocentric radius. However, 
a proper calculation of the tidal radius is still a challenge for dynamicists, including its 
definition itself. For our calculations we used the observed present day tidal radii (Table 1), 
which were obtained by fitting a single-mass King model to the observed cluster density 
profiles. 

Heating of stars, reversing the core collapse, is due to three-body binaries. They are 
included explicitly, without following their actual formation and evolution, according to the 
prescription by [Cohn (1985)| . Although [Ostriker (1985)| showed that the tidally captured 



binaries are probably more dynamically important for massive clusters than those formed 
through the tree-body interactions, the re-expansion phase following the core collapse is 
largely independent of the central heating source (Henon 1961; Goodman 1993). 



2.7. Philosophy of the Numerical Experiments 

Our goal in this paper is to investigate the importance of the different destruction 
processes on the overall evolution of the Galactic Globular Cluster System. We consider 
the evaporation process, and the disk and bulge gravitational shocks. AHO included also 
the dynamical friction in their orbit integrations, but we cannot model this effect at the 
moment. We rely on the AHO's results that the dynamical friction is not an important 
destruction mechanism for most clusters at the present time. 

Also, AHO evaluated the destruction rates for all of the mechanisms separately. Using 
the Fokker-Planck code we could investigate the effects only together, acting simultaneously. 
We hope this brings the numerical simulations closer to the reality, since in the real clusters 
all of these processes act coherently, thus "helping" each other. For example, the disk and 
bulge shock induced relaxation (eqs. [12 .2^) combines with the normal two-body relaxation 



and enhances the evaporation of stars through the tidal boundary. 



- 16 - 



Nevertheless, it is of prime interest to rank the processes in their role for the destruction 
of the clusters. In particular, AHO found that the bulge shock is more important mechanism 
than the disk shock. Also the newly discovered tidal shock relaxation (KO) has never been 
used in the detailed calculations of globular cluster evolution. How important are these 
processes? We try to answer this question using a "reduction" approach. We perform 
several sets of runs of the Fokker-Planck code, increasing the number of the processes 
allowed to act upon the cluster. The magnitude of each effect can then be estimated by 
subtracting the results of the run without that effect from the results of the run including 
the effect. Since we do not expect the effects to add in a linear fashion such an estimate will 
only approximately determine the strength of the effect. Thus we organized the following 
series: 0) evaporation only; evaporation + 1) disk shock, (AE) term only, no tidal shock 
relaxation; 2) disk shock, including the relaxation term {AE"^); 3) disk shock + bulge 
shock, without tidal shock relaxation; and finally 4) disk shock + bulge shock, including the 
relaxation. The last run includes all the processes we consider in this paper and represents 
the total destruction rate computed for the globular clusters. 

For each set of runs we repeat the calculations three times, allowing for different 
adiabatic conservation factors for the shock processes (sections p.5.2| and p.5.3|) . The first 



run includes the Spitzer correction, the second - Weinberg correction, and the last one 
assumes the impulse approximation with no adiabatic correction applied. 

Each run includes the Fokker-Planck simulations for all of the clusters in our sample, 
starting from the present time and ending with their total destruction. We use the observed 
concentrations and core radii of the clusters (Table 1) to model their current structure, 
which we approximate by a single-mass King model. The statistically assigned kinematic 
parameters {Rp, Ra, Vp] section |2.4| ) were used to estimate the amplitude of the tidal 
shocks. We followed evolution of the cluster up to a late stage near total destruction when 
the code breaks down due to numerical difficulties in recomputation the cluster potential. 
The remaining mass at that point is on average 8% (but no more than 11%) of the initial 
cluster mass. Destruction time td was extrapolated using the least-squares linear fit to 
the last 10 integration steps. This gives a more robust estimate of td than just linear 
extrapolation from the last couple points because many clusters suffer the gravothermal 
instability (e.g., [Goodman] 1993) at the late stages of their evolution. We return to this 



issue in the Results section. 

Given the destruction times for all the clusters in the sample, we obtain the destruction 
rate in units of a Hubble time 

^Hubble 10^° yr 
u = — — = — . (26) 

l^d T^d 

This definition agrees with the destruction rates of AHO, thus allowing a direct comparison. 
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Since the resulting distribution of the rates is broad and in general asymmetric around the 
mean, we choose to take a median z/ of the sample as a characteristic value, in accordance 
with AHO. Also, a few clusters from the sample have v ~ 10^ and therefore dominate in the 
mean and standard deviation. The standard error of the median was estimated as ([Kendall, 
Stuart, fc Ord 1987|) 

_ 1 

~ 2 /(/>)' ^^^^ 
where fiv) is the probability density distribution evaluated at the median point. For 
our sample we write ]{v) = AN/NAu, since / should be normalized to unity. Taking 
AA^ = O.bN^^"^, we define the two-sided errors 

af = iy{N/2 + N^/y2)-i>, (28a) 
err = u{N/2- N^/y2) -D, (28b) 

for the sorted sample of u. 



3. Results 

The results of our calculations for the individual clusters are presented in Table 3. 
Here we show the destruction rates for the two galactic models and the isotropic starting 
kinematics of the clusters (the anisotropic model is omitted for brevity). The Weinberg 
adiabatic correction is assumed throughout the Table. The first two columns identify the 
clusters and are the same as in Table 1. The third column gives the half- mass relaxation 
time in years. The remaining columns present the median destruction rates per Hubble 
time for each cluster in the sample. The first is the evaporation rate, followed by the runs 
including the tidal shocks: first and second order disk shock, and first order disk-|-bulge 
shocks for both the galactic models. The columns 7 and 10 give the total destruction rates 
corresponding to our present time globular clusters. The rest of this section discusses the 
statistical results for the sample. 



3.1. Evaporation 

We turn first to the evaporation of stars from clusters. In other words we perform a 
set of integrations where we allow only for ordinary two-body relaxation and losses over the 
tidal boundary. Although many studies have been devoted to normal two-body relaxation, 
we are unaware of any complete survey of a distribution of destruction times for different 
cluster parameters. Figure |^ shows the evaporation time in units of the initial relaxation 
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time versus cluster structural parameter, concentration c = log{Rt/ Rc), where Rt and Rc 
are the tidal and core radii, respectively. 

Almost all points on this Figure lie along one curve. This makes it very useful for 
quick estimates of the evaporation time given the initial parameters of a cluster. The 
least-squares fit gives 

^ = 1.290 X 10-1 - 1.170 X 10-1 c + 3.282 x 10-^ + '^•355x10-^ ^ 
Tev (c - 0.55)^ 

where Trh is the half-mass relaxation time, and Tev is the total evaporation time. Our fit is 
shown on Figure ^ as a solid line. 

This result suggests that the loosely bound clusters, with c < 0.65, are destroyed very 
fast (in units of the relaxation time). The most stable against relaxation are the clusters 
with c = 1.5 — 2, which survive for about 40 trh- But this number decreases to 20 trh as the 
concentration rises beyond 2. Similar behavior has been found by [Johnston^ (1993; his Fig. 
2)- 

The correlation between T^^ and Trh could be attributed to the fact that the King 
models (the initial condition for our calculations) belong a one-parameter family, and so can 
be described by the concentration parameter only. This holds only for the dimensionless 
quantities, such as the ratio Tev /Trh- The plot of the current relaxation time, expressed 
in years (Figure |^), does not show any obvious correlation with the concentration of the 
clusters. 



3.2. Gravitational Shocks 

Inclusion of the gravitational shocks speeds up the destruction of the clusters 
dramatically. Table ^ summarizes the results of our simulations when all the physical 
processes are acting on the clusters. The destruction rates are calculated as medians for the 
sample (see Section |2.7| ). In the OC galactic model the shocks almost double the destruction 
rate due to relaxation. Disk shocking ("old" first order) has little effect for both isotropic 
and anisotropic kinematic models, but the shock induced relaxation is more pronounced 
(0.09 in the median rate gain versus 0.017 in the former case). The case with the Spitzer 
correction is enhanced in the shock relaxation and bulge shocks because of the over-impulse 
increase of the adiabatic corrections (Section |2.5.2|) . The bulge shock is much stronger, and 
in the isotropic model it completely dominates over the disk shock. For the anisotropic 
model, the destruction associated with the bulge shock is reduced. Total destruction rate 
for the OC model is in the range 0.45 — 0.58. 
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In the BSS galactic model the disk is stronger, and there is a nuclear component that 
dominates the tidal shock over large range of radii (cf Figure Therefore the destruction 
rates are significantly increased by the shocks: they range from 0.63 to 0.86. 

Next, we investigate the dependence of the destruction rates on the current observed 
cluster position. Figure || shows the destruction rates associated with the several 
combinations of the physical mechanisms (Section |2.7| ) versus galactocentric radius, for the 
OC galactic model and isotropic kinematic model for the clusters (Section |2]^). We use 
here calculations with the Weinberg adiabatic correction. The left top panel, corresponding 
to the two-body relaxation, is largely determined by the selection of our sample, since 
relaxation is an internal process. However the steady increase in the rate with the decreasing 
distance to the galactic center is most likely explained by the growing tidal field that 
imposes tidal cutoff, and therefore removes stars from cluster faster. Going from the left 
top panel to the middle bottom we include more and more tidal shock processes. The disk 
in the OC model is relatively weak and does not change the distribution noticeably. Bulge 
shocks, on the other hand, enhance the destruction in the center considerably. Relaxation 
induced by the bulge shock is comparable to the first order effect. For clusters within 2 kpc 
of the center of the Galaxy, tidal shock relaxation dominates ordinary two-body relaxation 
in determining the cluster evolution. For such clusters core collapse occurs much faster and 
overall evolution proceeds on the shock time scale, when it is shorter than the relaxation 
time. During the cluster contraction, the relative importance of the tidal shock relaxation 
to ordinary two-body relaxation decreases and the final stage of core collapse is described 
by the self-similar solutions of |Henon 1961i After core collapse the cluster is still so highly 



concentrated that the density profile is close to an isothermal sphere. Slow expansion along 
with the tidal stripping of stars leads to final dissolution of the cluster. 

A similar plot for the anisotropic kinematic model is given in Figure ^. Figures ^ 
and |ri] show the distribution for the BSS galactic model and the isotropic and anisotropic 
kinematics, respectively. 

To emphasize the relative importance of each of the destruction mechanisms, we define 
the differential rates by subtracting the destruction rates obtained without that process 
from the run including the process. Namely, first-order disk shock ("diskl") is "evaporation 
-|- diskl" - "evaporation" , second order disk shock relaxation term ( "disk2" ) is "evaporation 
+ disk" - "evaporation -|- diskl", first order bulge shock ("bulgel") is "evaporation -|- 
diskl -|- bulgel" - "evaporation -|- diskl", and the bulge relaxation is "evaporation -|- 
disk -|- bulge" - "evaporation -|- disk" - "bulgel" . Differential rates for the two Galactic 



and kinematic models are presented in Figures - |15| For the reasons noted earlier 



(nonlinearity), these figures must be considered indicative but not exact. 
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We present the histograms of the distributions of the destruction rates in Figures 16 



T9|. Arrows at the top show median of the distribution. 



We should however be cautious using these differential results, since all those 
mechanisms act together and the final result is not a direct sum of the single processes. 
For example, relaxation is greatly enhanced by the tidal shock relaxation and core collapse 
occurs faster (in some cases clusters collapsed even when the ordinary relaxation was not 
enough to drive the contraction). 



3.3. Vital Diagram for Globular Clusters 

We do not know what was the initial distribution of globular clusters in our Galaxy. 
But we can imagine that they occupied some volume in a given parameter space. All the 
physical mechanisms considered above tend to destroy clusters with time, and thus to 
reduce the allowed volume. They superimpose particular boundaries which distinguish the 
present day clusters from those being already dissolved (or never formed). [I'all fc Rees| 
(1977)1 considered the cluster mass versus their typical size (half-mass radius Rh) diagram. 



The authors included evaporation, disk shock and tidal shock heating (for the latter they 
assume the tidal interaction with the clusters themselves rather than with the bulge). 
These processes cut a triangle on the Rh — M plane, containing the observed clusters pretty 
well. Ostriker (1975, unpublished) and [Caputo fc Castellani (1984] also included dynamical 



friction, which excludes the very massive clusters. They noted also that the strength of the 
destruction mechanisms vary with galactocentric radius, and therefore the allowed space 
depends on the position of a cluster. 

We constructed such a vital diagram for the Galactic GCS using our sample and results 
of the computations of its evolution. All of the processes used in the simulations, as well as 
the dynamical friction, participate in the diagram. The vital boundary is defined in such a 
way, that the sum of all the destruction rates is equal to the inverse Hubble time: 

1 _ 1 ^ 1 ^ 1 (^3Q^ 

'^Hubble ^ev ^sh ^df 

where t^v, tsh, and t^f are the time-scales over which a cluster would be destroyed 
by the given process alone, for the evaporation, disk and bulge shock combined, and 
dynamical friction, respectively. Note that throughout this paper we adopted for simplicity 
^Hubble = 10^° yr. 



According to the calculations reported in Section 

tev ~ 30trh, (31) 
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where trh is the half-mass relaxation time (eq. 0). We plan to perform a similar 
comprehensive analysis for the gravitational shocks. At present we use 

tsh = (32) 

^'disk + i^bulge 

Thus the shock boundary should be considered as a stronger limit. 



We have not considered the effects of dynamical friction in details in this paper. Pinney 
fc Tremaine (198"7]| estimated the time for a cluster to lose its momentum and fall to the 
Galactic center (their eq. 7-26): 



2.64 X 10"yr f R V f V, \ f Mq \ 
InA 1^2kpcy' l250kms-V M j ' 

where InA is the Coulomb logarithm, R is the initial galactocentric distance, Vc is the 
cluster circular speed, and M is its mass. 



The cluster vital diagram for the OC galactic model is shown in Figures 20, 21 



The disgram for the BSS model is in Figure 22, 23 



One can see that a significant fraction of the clusters lies outside of the "surviving" 
boundary. This supports our conclusion that many clusters will be destroyed within the 
next Hubble time. 



4. Discussion 

Using the Fokker-Planck simulations and the sample of the globular clusters we 
estimated a current destruction rate per Hubble time. These results are applicable for the 
present day clusters evolving forward in time for the next Hubble time. On the other hand, 
it is of prime interest to extract any possible information regarding the past evolution of the 
clusters from the time of their formation up to now. We try to construct a simple model 
for the initial distribution of the clusters and to test it against the current destruction rate. 
We try then to answer the question raised in the Introduction: how many of the globular 
clusters may have been destroyed in the history of our Galaxy. 

We will not attempt to propose a mechanism for the formation of globular clusters. 
An example of a possible formation scenario is given by [Fall fc Rees (1985]] . Instead, we 



assume a life-time function for the globular clusters, free of any particular assumptions. 
Let td{t) be a time remaining to the total destruction of a given cluster at epoch t. We 
then define /(t^; t) dt^ to be the number of clusters with the destruction time in the interval 
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\td-, td + dtd] at time t. Thus, if we treat almost all clusters as made in a short time interval 
at the formation of the Galaxy (an obviously gross over-simplification), then f{td',0) dtd 
gives the initial distribution of cluster lifetimes. The normalization is 

N{t) = / fiU;t)dU, (34) 





where N{t) is the number of clusters at epoch t. We would like to advance our function / in 
time starting from their formation, again without a detailed prescription for the evolution 
of the cluster sample. We simply assume that all the clusters were formed at the same time, 
and that this time to ~ is very small compare to the present Hubble time. Then number 
of clusters surviving at the time t is just the number of initial clusters with td > t. Thus we 
have the following relation: 

f{td;t) = f{td + t;to), (35) 

where we have neglected the contribution of to ^ ^ iii the first argument of the function /. 
We call fitd]to) = fiitd) the initial distribution of the globular clusters. Integrated over all 
destruction times td it gives the initial number of globular clusters N{t = to) = Ni formed 
in our Galaxy. 

We define also mean td and median t^ destruction times according to 

r (.\ _ /o°° id fjtd; t) dtd _ /o°° td f{td\ t) dtd , , 

k^f{td,t)dtd - Nit) ' ^'^^ 

1 ^ /o*'" f{td]t)dtd ^ /q- f{td]t)dtd 

2 Jo°° f{td;t)dtd N{t) • ^ ^ 

Now we choose two functional forms for fi{td), which we test against the distribution 
of the destruction rates obtained in Section |^. The simplest one is that with a constant 
mean destruction time for all clusters. It assumes an exponential form 

Mtd) = Cie-"*-^. (38) 

It is easy to show that the mean and median for this distribution are 

idi = -, (39) 
a 

ln2 , , 

tml = • 40 

a 

Thus the mean destruction rate does not depend on time, and the evaporation of the 
clusters in this case is similar to radioactive decay. Given a current vital rate we can always 
calculate how many clusters survived from the beginning. N{t) goes exponentially to zero 

NUt) = —e'''\ (41) 
a 
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The other function / we consider is a scale free power-law 

UU) = C^t-/. (42) 

This distribution is not strictly normalizable since the total number of clusters iVj diverges 
as to approaches zero 

iV2(t) = ^^^-^ g>l. (43) 

q-l 

However we can successfully apply it to the present time. The corresponding mean and 
median are 

id2{t) = q>2, (44) 

q-2 

tm2it) = t(2^-l). (45) 

Note that both these quantities are proportional to the time of observation. Thus it seems 
impossible to determine the initial cluster population for this distribution given the current 
rate. Fortunately, we can compare the shape of the distribution with the current profile 
to choose between the two distributions. In addition we note a very important difference 
consequent to the two hypothesized forms for the initial distribution. In the second (power 
law) case we should expect that id ~ t; whenever we look, the time to destruction for the 
clusters that remain is of the same order as the age of the existing sample. In the first case 
this might occur, but it would require a coincidence between the initial typical time to 
destruction and the later point in time when an observer examines the system. The most 
important point to be made in this paper is that the numbers in the last line of Table are 
of order unity. This is consistent with the power law assumption and allows the possibility 
of large fractional destruction. 



On Figure we plot a histogram of the destruction rates for our sample of 119 clusters 



for the OC galactic model. The dotted region gives the number of clusters per logarithmic 
bin in t^- Two panels correspond to the results obtained for the isotropic kinematic model 
with the Spitzer and Weinberg adiabatic corrections, respectively. We use the medians 
for the two cases to determine the shape of the two models for the destruction rates. The 
dashed curve corresponds to the exponential model, and the solid line is for the power-law 
one. Both lines are normalized to the present number of clusters, iV(tHubbie)- 

Figure shows the distribution for the BSS galactic model. 

We see from this plot that both models are consistent with the actual distribution. 
However, in our view, the power-law case better represents the overall extended shape of 
the histogram. The exponential model is more concentrated around the median value. 
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and falls off very rapidly for large t^. It predicts also much larger number of clusters in 
the center of the diagram than follows from the Fokker-Planck calculations. Therefore we 
favor the second choice, namely, the power-law distribution of the destruction times. Based 
on the results obtained for the OC galactic model and Weinberg adiabatic correction, 
power-law index is g ~ 1.59 — 1.64 (two limits are for the anisotropic and isotropic kinematic 
distributions, respectively). For the BSS model we find q ~ 1.73 — 1.82. Thus we conclude 
that the power-law model with g = 1.6 — 1.8 gives naturally the present distribution of the 
globular cluster destruction times. 

From the last sentence immediately follows that the initial cluster population could 
have been much larger than the present one. The present day characteristics do not provide 
a definite number, though they do prefer this scenario. Thus it is quite possible that a 
very significant fraction of the original globular clusters have been already destroyed. Their 
remnants might constitute now the inner spherical stellar component of our Galaxy, i.e. 
the bulge. We note also from Figures | and |l^ that clusters close to the Galactic centers 
are most heavily attacked by the bulge shock, and therefore most of the dissolved clusters 
could be in the very central region of the Galaxy. A strong bulge component present in 
the Ostriker & Caldwell galactic model could be enhanced by the input from the destroyed 
clusters. Note, that we do not require the existence of the bulge from the first days of the 
Galactic history in order to destroy the clusters. As follows from the power-law distribution, 
most of the initially formed clusters had relatively short life-times, and presumably were 
not very massive. Thus the main mechanism driving their dissolution could have been 
internal relaxation. Though, if present in some form, the central Galactic component would 
be very efficient in the cluster destruction. In any case the overall potential of the Galaxy 
is changed only very slightly if we were to imagine that all bulge and spheroid stars were 
put back into globular clusters with appropriate orbits. Thus even the bulge shock tidal 
evaporation processes would be essentially unchanged, if we were to put all stars in the 
quasispherical distribution back into clusters. 

Observational data on stellar populations indicate that the spheroid of the Galaxy 
is kinematically and chemically distinct from the disk ( [Norris fc Ryan 1991|) . Halo stars 
are old (with an age slightly exceeding our adopted Hubble time of 10^° yr) and metal 
poor. Because the density of the protogalactic gas cloud is not high enough to form 
stars efficiently, the oldest stars in the Galaxy are likely to have been formed in clusters 
perhaps not too different from the present population of globular clusters. Now we see 
most of the metal poor stars in the halo, which suggest that most of them could be left 
over from the disrupted clusters. [Harris (1991)| proposed two arguments against this idea: 



1) orbits of the present globular clusters are more isotropically distributed than those of 
spheroid stars; 2) the clusters have systematically lower metallicity than the field stars, and 
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become more metal poor with increasing galactocentric distance R. But AHO and |Lee fc 



Goodman] (1995) pointed out a possible solution for the former problem. The first of these 
arguments can be moderated if we consider the initial cluster distribution. AHO argued 
that the kinematic model better describing the current population is the isotropic one (see 
Section |2.3| ) with the orbits getting more and more radial as R increases. The destruction 
mechanisms are the strongest in the central few kiloparsecs, especially the bulge shock. 
Thus the clusters in the most elongated orbits with small perigalactic distance Rp would 
be destroyed first. Their remnants would populate the spheroid as wee see it now, with 
the preferentially radial orbits of stars. Surviving clusters, on the other hand, have more 
isotropic orbits consistent with the observations. 

To resolve the problem of metallicities, we refer to |van den Bergh (1995)| , who found 
that [Fe/H] correlates somewhat more strongly with perigalacticon Rp than with their 
current position R. Thus the metal rich clusters with smaller Rp, which would be destroyed 
first, would enrich the halo stellar population. Another solution comes from the observed 
age difference in the cluster population. The formation period could have been extended up 
to 2 — 5 Gyr (see references in |N orris fc Ryan 1991|) . Thus the very first low mass clusters 
could dissolve due to the internal relaxation, but their stars would enrich the primordial 
gas for the next generation of clusters. Our power-law hypothesis for the distribution of the 
clusters predicts that the mean (and median) destruction time is proportional to the time 
elapsed since the formation, so that newly formed clusters would be again most susceptible 
to disruption. If they were destroyed within the remaining 7—10 Gyr, their stars are the 
population II halo stars that we see in our Galaxy. 



5. Conclusions 

We have used the Fokker-Planck code to investigate the destruction rate of globular 
clusters in our Galaxy. We applied two forms of the adiabatic correction for gravitational 
shocks and found that the median results do not depend much of the particular form of 
the correction. The current destruction rate for the sample is about 0.5 — 0.9 per 10^° 
yr (depending on the Galactic and kinematic models), which implies that more than half 
of the "present clusters is to be destroyed within the next Hubble time. This estimate is 
approximately a factor of ten higher than that obtained by AHO. There are two principal 
reasons for the change. First, our Fokker-Planck detailed calculations for each cluster give 
systematically larger rates of two body relaxation, and core collapse than did the essentially 
time scale arguments of AHO. Secondly, the new tidal shock relaxation process described 
by KO further reduces lifetimes by a significant amount. 
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Trying to understand the original population of Galactic globular clusters, we 
considered two possible models for the distribution of the cluster lifetimes. Both of them 
can be normalized to the median present destruction rate. We favor the power-law model 
on basis of a shape of the rate distribution as it naturally explains the fact that the current 
median time to destruction is comparable to the present mean cluster age and also because 
the predicted distribution of cluster destruction times provides a reasonable match to 
observations. The power-law distribution allows a much larger number of the clusters to 
have been formed initially than is currently observed, and allows the possibility that the 
debris of the early disrupted clusters might have formed the much of spheroid of our Galaxy. 

Surdin [1995^ has investigated the possibility of populating the stellar halo by remnants 
of the destroyed star cluster (open and globular). He comes to the conclusion that much 
larger number of low-mass (10^ — IO^Mq) clusters than observed now is required to match 
the mass of the spheroid. This conforms to our result, since low mass (and hence, weakly 
concentrated) clusters have short lifetime. All those non-observed clusters have to be 
destroyed before the present time. 

Finally we note that the inclusion of the mass spectrum in the Fokker-Planck models 
would strongly enhance the relaxation and core collapse, and ultimately speed up the 
dissolution of the clusters. This could also increase the destruction rate by a significant 
amount. 

The destruction rates for the individual clusters are available electronically upon 
request. 
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FIGURE CAPTIONS 

Fig. 1. — Histogram of the orbit eccentricity distribution in the OC galactic model, for the 
isotropic and anisotropic kinematic models. 

Fig. 2. — Histogram of the orbit eccentricity distribution in the BSS galactic model, for the 
isotropic and anisotropic kinematic models. 

Fig. 3. — Comparison of the Spitzer [eq. (|10D] and Weinberg [eq. (|TlD ] adiabatic corrections 
for the disk shock, i. e., reductions in the energy change due to the shock because of the 
conservation of adiabatic invariants of stellar orbits inside the clusters, versus adiabatic 
parameter [eq. (|)]. 

Fig. 4. — Correction factors versus perigalactic distance for the bulge shock energy change 
due to the extended mass distribution for the bulge. Our correction [eq. (p!8D] is shown by 
solid curve, and that from Aguilar & White (1985), and subsequently AHO, is given by dots. 
Triangles show the normalized bulge mass profile for comparison. If the bulge were treated 
as a point mass, that would have been the factor. 

Fig. 5. — Same as Figure |], but for the bulge shock. The adiabatic parameter is given by 
equation (|2T|). 

Fig. 6. — Evaporation time in units of initial relaxation time versus cluster concentration. 
Dots are results of the Fokker-Planck calculations with the observational data as input. Solid 
line - our fit, equation (P^). 

Fig. 7. — Relaxation time at the half-mass radius versus cluster concentration for our sample 
of globular clusters. 

Fig. 8. — Mean probability per Hubble time for cluster destruction due to evaporation, 
and the evaporation plus gravitational shocks: disk shock without relaxation, disk shock 
including relaxation term, disk+bulge shocks without relaxation, and disk+bulge shock 
including relaxation, as a function of present day galactocentric distance. The error bars 
indicate the interquartile range of the distribution. Weinberg adiabatic corrections are used. 
Galactic model is OC, and kinematic model is isotropic. 

Fig. 9. — Same as Figure |^, but for the anisotropic kinematic model. 

Fig. 10. — Same as Figure ||, but for the BSS galactic model. 

Fig. 11. — Same as Figure |^, but for the BSS galactic model. 
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Fig. 12. — Differential destruction rate for different processes: evaporation, disk sliock witli 
and witliout induced relaxation ("diskl" and "disk", respectively), and bulge sliock with and 
without relaxation term ("bulgel" and "bulge"). See text for the definition of the differential 
rate. The error bars indicate the interquartile range of the distribution. Weinberg adiabatic 
corrections are used. Galactic model is OC, and kinematic model is isotropic. 



Fig. 13. — Same as Figure 12, but for the anisotropic kinematic model. 



Fig. 14. — Same as Figure 12, but for the BSS galactic model. 



Fig. 15. — Same as Figure |T3|, but for the BSS galactic model. 

Fig. 16. — Histogram distribution of the destruction rates from Figure ^j. Weinberg adiabatic 
correction is used. Galactic model is OC, and kinematic model is isotropic. 

Fig. 17. — Histogram distribution of the destruction rates from Figure Weinberg adiabatic 
correction is used. Galactic model is OC, and kinematic model is anisotropic. 

Fig. 18. — Histogram distribution of the destruction rates from Figure 0. Weinberg 
adiabatic correction is used. Galactic model is BSS, and kinematic model is isotropic. 

Fig. 19. — Histogram distribution of the destruction rates from Figure 0. Weinberg 
adiabatic correction is used. Galactic model is BSS, and kinematic model is anisotropic. 



Fig. 20. — Vital diagram for the Galactic globular clusters. Mass-radius plane is restricted 
by three destructing processes: relaxation, tidal shocks, and dynamical friction. Galactic 
model is OC, and the kinematic model is isotropic. 



Fig. 21. — Same as Figure |20|, but for the anisotropic kinematic model. 



Fig. 22. — Same as Figure |20|, but for the BSS galactic model. 



Fig. 23. — Same as Figure |21], but for the BSS galactic model. 

Fig. 24. — Distribution of the destruction times for our sample for the OC galactic model. 
All destruction processes are included. Solid line - power-law distribution, normalized to the 
present time, dashes - normalized exponential model. 



Fig. 25. — Same as Figure 24, but for the BSS galactic model. 
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Table 2. Kinematic parameters of Galactic Globular Cluster System 



Reference 


Vrot 


0"los 




(km s~^) 


(km s~^) 


Frenk & White (1980) 


60 ±26 


118 ±20^ 


Thomas (1989) 


65 ±18 


110 ±7 


This paper 


60 ±21 


119 ±39 



'^The uncertainty of the velocity dispersion is taken as an average 
over radial bins of Table 1 from [White fc Frenk (1983) . 
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Table 4. Median Destruction Rates 







Isotropic- 


Anisotropic 


Destruction Process 


Spitzer 


Weinberg No correction 


Spitzer Weinberg 



Ostriker & Caldwell model 



Evaporation 0.290 1°-°?° 0.290 1°-°?° 0.290 1°-°?^ 0.290 1°-°?° 0.290 l°-°?3 

+ Disk 1st 0.307 ;°:°!? 0.311 ;°:°^^ 0.312 i°:°i? 0.307 ;°:°li 0.309;°;°?^ 

+ Disk 0.397 1°:°!^ 0.391 in^g 0.394 1°:°!^ 0.393 1^?! 0.365 1^1 

+ Disk 1st + Bulge 1st 0.467 l°:°i 0.455 1°;°^^ 0.601 t°:°^? 0.398 t°:°i 0.397 1°:°?? 

+ Disk + Bulge 0.583 l°:Ji 0.518 1°:J°« 0.905 1°:°^1 0.490 l°:Ji? 0.453 1°:J^^ 



Bahcall, Schmidt & Soneira model 



Evaporation 0.290 +°-°?3 0.290 +°-°?3 0.290 +°-°?3 0.290 +°-°?3 0.290 +°-°?° 

+ Disk 1st 0.495 l°:°i 0.498 1°:°?^ 0.502 1°;^^? 0.395 l^fr 0.395 1°;°^^ 

+ Disk 0.654 ;°:°SI 0.660 i°:°i 0.818 1°:°^^ 0.511 ;°:J1^ 0.493;°:°^^ 

+ Disk 1st + Bulge 1st 0.627 toltt 0.628 Ig-}!^ 1.026 0.520 1^54 0.480 l^lg 

+ Disk + Bulge 0.863 ;°;}i 0.752 1°:?!^ 1.786 1°;^^^ 0.712 1°;^°? 0.633 l°:°i 



Note. — Rates are in fraction of sample per Hubble time (10^° yr). Total number of clusters in 
sample: 119. 
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TABLE 1 
Globular Cluster Parameters 



ID 



Name 



D 



R 



Rn 



Rt 



M 



Vr 



Rcf 



> 

O 
m 
o 

ON 
Oh 

6 



















foe) 


foe) 




(km s 




1 


NGC0104 


305.90 


-44.89 


4.6 


7.8 


2 


04 


0.50 


54.8 


1.45e+06 


-19.0 


1 


2 


NGC0288 


152.28 


-89.38 


8.3 


11.9 





96 


3.47 


31.6 


l.llc+05 


-46.0 


1 


3 


NGC0362 


301.53 


-46.25 


8.5 


9.6 


1 


94 


0.43 


37.5 


3.78e-h05 


223.2 


1 


4 


NGC1261 


270.54 


-52.13 


15.9 


18.0 


1 


27 


1.82 


33.9 


3.26e-h05 


51.0 


9 


5 


Pall 


130.07 


19.03 


13.7 


20.1 


1 


50 


0.60 


19.0 


2.54C+03 


3.0 


9 


6 


AMI 


258.36 


-48.47 


124.7 


126.1 


1 


23 


4.57 


77.6 


1.81e-h04 


116.0 


11 


7 


Eridanus 


218.11 


-41.33 


79.5 


84.8 


1 


10 


5.89 


74.2 


2.30e-h04 


-21.0 


10 


8 


Pal2 


170.53 


-9.07 


13.6 


21.9 


1 


40 


0.60 


15.1 


1.74e-H05 


-133.0 


9 


9 


NGC1851 


244.51 


-35.04 


12.3 


17.3 


2 


24 


0.28 


48.7 


5.61e-h05 


320.5 


1 


10 


NGC1904 


227.23 


-29.35 


12.9 


19.2 


1 


72 


0.60 


31.5 


3.57e-h05 


203.4 


1 


11 


NGC2298 


245.63 


-16.01 


10.0 


15.5 


1 


28 


1.00 


19.1 


7.20e-t-04 


150.4 


3 


12 


NGC2419 


180.37 


25.24 


83.2 


91.0 


1 


40 


8.51 


213.8 


1.60e-h06 


-20.0 


1 


13 


NGC2808 


282.19 


-11.25 


9.1 


11.1 


1 


77 


0.71 


41.8 


1.32e-H06 


98.0 


9 


14 


Pal3 


240.14 


41.86 


90.3 


93.8 


1 


00 


12.59 


125.9 


6.38e+04 


84.0 


9 


15 


NGC3201 


277.23 


8.64 


5.1 


9.3 


1 


31 


2.09 


42.7 


1.95e-h05 


494.6 


1 


16 


Pal4 


202.31 


71.80 


98.2 


101.0 





78 


15.85 


95.5 


5.41e-h04 


75.0 


9 


17 


NGC4147 


252.85 


77.19 


18.7 


21.0 


1 


80 


0.55 


34.7 


6.69e-h04 


183.0 


1 


18 


NGC4372 


300.99 


-9.88 


5.2 


7.4 


1 


30 


2.63 


52.5 


3.20e-|-05 


49.0 


9 


19 


RuplOG 


300.89 


11.67 


19.7 


17.1 





70 


5.75 


28.8 


8.42e-f04 


-44.0 


7 


20 


NGC4590 


299.63 


36.05 


9.3 


9.8 


1 


64 


1.86 


81.2 


3.06e-|-05 


-95.1 


1 


21 


NGC4833 


303.61 


-8.01 


5.8 


7.2 


1 


25 


1.70 


30.2 


9.49e-f04 


194.0 


9 


22 


NGC5024 


332.96 


79.77 


18.5 


19.1 


1 


78 


2.00 


120.5 


7.33e-t-05 


-80.0 


9 


23 


NGC5053 


335.69 


78.94 


16.1 


16.9 





82 


10.47 


69.2 


1.66e-|-05 


42.8 
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24 
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309.10 


14.97 
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6.8 
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232.2 
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25 
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42.21 
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10.2 
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43 
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5.5 


4.7 


1 


38 


1.07 


25.7 


4.93e-h05 


-42.6 


1 


47 


NGC6229 


73.64 


40.31 


31.6 


30.9 


1 


61 


1.23 


50.1 


4.46e-H05 


-154.2 


8 


48 


NGC6235 


358.92 


13.52 


9.6 


2.4 


1 


33 


1.00 


21.4 


1.71e-f05 


86.0 


9 



TABLE 1— Continued 



ID 


"XT 

Name 


7 

(deg) 


7 



(deg) 


D 
(kpc) 


R 
(kpc) 


c 


Kc 

(pc) 


Rt 

(pc) 


M 


T r 

Vr 

(km s ) 


Ret 


49 


NOC6254 


15.14 


23.08 


4.3 


5.1 


1.40 


1.07 


26.9 


2 25P-I-05 


75.5 


1 


50 


Nr;C6256 


347.79 


3.31 


10.2 


2.7 


2.50 


0.07 


22.1 




-104.7 


1 


51 


Pall5 


18.87 


24.30 


41.3 


34.2 


0.60 


14.45 


57.5 




68.9 


1 


52 


NGC6266 


353.57 


7.32 


5.6 


3.1 


1.70 


0.29 


14.5 




-71.9 


1 


53 




356.87 


9.38 


10.9 


2.9 


1.53 


1.35 


45.7 


-L .(JVJL. 1 \j\J 


131.0 


9 


54 


NOC6284 


358 35 


9.94 


11.8 


3.8 


2.50 


0.25 


79.1 


2 1 7P-I-05 


27.4 


1 


55 


NOC6287 


0.13 


11.02 


6.7 


2.3 


1.60 


0.52 


20.7 


-L • v/OC 1 \)^J 


-208.0 


9 


56 


i 1 \J \J ^ KJ 


357.62 


7.83 


8.4 


1.2 


2.50 


0.11 


34.8 


2 34o+05 


-148.0 


1 


57 


NC1C6304 


355 S3 


5.38 


6.0 


2.7 


1.80 


0.36 


22.7 


1 qsp-i-n5 


-98.0 


8 


58 


NGC6316 


357.18 


5.76 


12.8 


4.5 


1.55 


0.62 


22.0 


1 n2e-l-06 

-L • \J s^yD 1 \J\J 


76.0 


9 


59 


NGC6325 


0.97 


8.00 


6.0 


2.7 


2.50 


0.05 


15.8 


9 57o+04 


30.9 


1 


60 


NGC6341 


68.34 


34.86 


7.5 


9.5 


1.81 


0.51 


32.9 


3 fi4p_|_n5 


-120.5 


8 


61 




5.54 


10.71 


7.4 


2.0 


1.15 


1.26 


17.8 


2 89P-I-05 


224.7 


8 


62 


NnCfi342 


4.90 


9.72 


12.1 


4.1 


2.50 


0.16 


50.6 


1 95o+05 


117.9 


1 


63 


NCC6356 


6.72 


10.22 


16.6 


8.5 


1.54 


1.12 


38.8 


« lOpinc; 


31.6 


8 


64 


NGC6355 


359 58 


5.43 


7.0 


1.7 


2.50 


0.10 


31.6 


3 fi7p-|.n5 


-184.0 


9 


65 


NGC6352 


341.42 


-7.17 


6.1 


3.5 


1.10 


1.48 


18.6 


1 22P-I-05 


-118.0 


9 


66 


NGC6366 


18.41 


16.04 


4.0 


5.1 


0.92 


2.14 


17.8 


■? q2p-i-n4 


-122.6 


1 


67 


NGC6362 


325 56 


-17.57 


7.8 


5.4 


1.10 


2.95 


37.1 


-L • J- 1 1 \)^J 


-13.3 


1 


68 


NGC6388 


345.56 


-6.74 


11.0 


3.7 


1.70 


0.40 


20.0 


1 50o+0fi 

-L .rJ\JKj \ \j\J 


77.0 


9 


69 


NGC6402 


21.32 


14.80 


10.1 


4.4 


1.60 


2.45 


97.5 


1 23e-l-06 

X • £j%J\Zi 1 \J\J 


-123.4 


8 


70 


NGC6401 


3.45 


3.98 


6.3 


2.3 


1.69 


0.45 


22.0 


1 23e-l-06 

J- • S^KtKj 1 \J\J 


-62.0 


9 


71 


NGC6397 

J. 1 v_J V_.' \J t/ 1 


338.17 


-11.96 


2.2 


6.5 


2.50 


0.03 


9.5 


1 5qp_|_n5 


17.7 


1 


72 


NGC6426 


28.09 


16.23 


17.5 


11.3 


1.70 


1.35 


67.7 


q 4qp_i_04 


-162.0 


9 


73 


NGC6440 


7.73 


3.80 


7.1 


1.8 


1.70 


0.26 


13.0 


5 72p-|-05 


-83.0 


9 


74 


NGC6441 


353 53 


-5.00 


10.7 


2.5 


1.85 


0.35 


24.8 


1 30o+0fi 


16.0 


1 


75 


NGC6453 


355.72 


-3.87 


10.1 


1.8 


2.50 


0.20 


63.2 


1 ,S2e-l-05 


-84.0 


9 


76 


NGG6496 


348.02 


-10.01 


5.8 


3.3 


0.70 


1.74 


8.7 


4 63p+f)4 


-95.0 


9 


77 


NGC651 7 


19.23 


6.76 


6.1 


3.5 


1.82 


0.11 


7.3 


1 82P-I-05 


-37.0 


9 


78 


NGC6522 


1.02 


-3.93 


7.2 


1.4 


2.50 


0.11 


34.8 


5 q3p_l_n4 


-10.4 


1 


79 




27.18 


10.44 


7.0 


4.1 


1.30 


0.83 


16.6 


5 q3p-i-n4 


-215.3 


1 


80 


NGC6528 


1.14 


-4.17 


6.6 


2.0 


2.29 


0.17 


33.1 


q 3ip_|_n4 


160.0 


9 


81 


NGC6539 


20.80 


6.78 


4.0 


5.0 


1.60 


0.63 


25.1 


1 fiOp-l-05 


-52.0 


9 


82 


NGC6544 


5.84 


-2.20 


2.6 


5.9 


2.50 


0.03 


9.5 


1 8f)p+n5 


-12.0 


8 


83 


NGC6541 


349.29 


-11.18 


6.7 


2.7 


2.00 


0.58 


58.0 


4 67p+05 


-152.8 


8 


84 


NGC6553 


5.25 


-3.02 


3.6 


4.9 


1.17 


0.56 


8.3 




-24.0 


9 


85 


NGC6558 


0.20 


-6.03 


8.9 


1.0 


2.50 


0.09 


28.5 


2 41p+n5 


-198.9 


1 


86 


NGC6569 


0.48 


-6.68 


8.8 


1.0 


1.27 


0.95 


17.7 


3 85p+05 


-26.0 


9 


87 


NGC6584 


342.14 


-16.41 


9.5 


3.9 


1.20 


1.66 


26.3 


2 iqp-i-05 


180.0 


8 


88 




2 7Q 


-7 Q1 

i .t/ J- 


8 1 

O. -L 


1 3 


2 50 


15 


47 4 


3 32p-I-05 


54.2 


\ 


89 


NGC6626 


7.80 


-5.58 


6.0 


2.8 


1.67 


0.42 


19.6 


4.42C+05 


15.5 


1 


90 


NGC6638 


7.90 


-7.15 


8.4 


1.6 


1.40 


0.65 


16.3 


1.12e-F05 


-14.0 


8 


91 


NGC6637 


1.72 


-10.27 


10.2 


2.4 


1.39 


1.00 


24.5 


3.44e-F05 


50.1 


8 


92 


NGC6642 


9.81 


-6.44 


7.9 


1.8 


1.99 


0.24 


23.5 


1.26e+05 


-47.0 


9 


93 


NGC6652 


1.53 


-11.38 


14.3 


6.2 


1.80 


0.30 


18.9 


2.59e-F05 


-124.2 


8 


94 


NGC6656 


9.89 


-7.55 


3.0 


5.6 


1.31 


1.23 


25.1 


5.36e-F05 


-148.8 


13 


95 


Pal8 


14.11 


-6.80 


28.6 


20.6 


0.75 


3.31 


18.6 


2.21e-h05 


-38.0 


15 


96 


NGC6681 


2.85 


-12.51 


9.2 


2.1 


2.50 


0.07 


22.1 


1.89e-F05 


221.1 


1 



TABLE 1— Continued 



ID 


Name 


I 


b 


D 


R 


c 


Rc 


Rt 


M 


Vr 


Ref 






(deg) 


(deg) 


(kpc) 


(kpc) 




(pc) 


(pc) 




(km s-i) 




97 


NGC6712 


25.35 


-4.32 


6.8 


3.8 


0.90 


1.86 


14.8 


2.45e-l-05 


-107.5 


1 


98 


NGC6715 


5.61 


-14.09 


21.1 


13.1 


1.84 


0.66 


45.7 


1.45e-|-06 


142.0 


3 


99 


NGC6717 


12.88 


-10.90 


7.3 


2.6 


2.07 


0.18 


21.1 


1.23e-|-05 


-6.0 


9 


100 


NGC6723 


0.07 


-17.30 


8.6 


2.6 


1.05 


2.40 


26.9 


3.74e-|-05 


-79.0 


9 


101 


NGC6752 


336.50 


-25.63 


4.2 


5.5 


2.50 


0.21 


66.4 


3.64e-|-05 


-32.1 


1 


102 


NGC6760 


36.11 


-3.92 


6.1 


5.1 


1.59 


0.59 


23.0 


1.38e-|-05 


-28.0 


9 


103 


Terzan7 


3.39 


-20.07 


36.7 


28.9 


1.08 


6.46 


77.7 


6.15e-h04 


162.0 


2 


104 


NGC6779 


62.66 


8.34 


9.7 


9.6 


1.37 


1.05 


24.6 


1.89e-f05 


-135.9 


1 


105 


Arp2 


8.55 


-20.79 


28.2 


20.6 


0.90 


13.18 


104.7 


3.26e-|-04 


119.0 


2 


106 


NGC6809 


8.79 


-23.27 


4.8 


4.6 


0.76 


3.98 


22.9 


2.41e-t-05 


175.5 


1 


107 


Terzan8 


5.76 


-24.56 


23.0 


15.7 


0.60 


6.61 


26.3 


2.50e-h04 


130.0 


2 


108 


Palll 


31.80 


-15.58 


13.2 


7.9 


0.69 


7.76 


38.0 


1.38e-|-05 


-68.0 


8 


109 


NGC6838 


56.74 


-4.56 


3.9 


7.2 


1.15 


0.72 


10.2 


3.67e-|-04 


-23.0 


1 


110 


NGC6864 


20.30 


-25.75 


18.7 


12.4 


1.88 


0.51 


38.7 


4.89e-F05 


-195.2 


8 


111 


NGC6934 


52.10 


-18.89 


14.8 


12.0 


1.53 


1.07 


36.3 


2.15e-h05 


-412.2 


1 


112 


NGC6981 


35.16 


-32.68 


18.1 


13.7 


1.23 


2.82 


47.9 


1.81e+05 


-309.0 


9 


113 


NGC7006 


63.77 


-19.41 


38.8 


36.1 


1.42 


2.75 


72.3 


2.52e-F05 


-364.0 


9 


114 


NGC7078 


65.01 


-27.31 


10.4 


10.7 


2.50 


0.21 


66.4 


9.84e-|-05 


-107.1 


12 


115 


NGC7089 


53.38 


-35.78 


11.8 


10.7 


1.80 


1.17 


73.8 


8.81C+05 


-3.1 


1 


116 


NGC7099 


27.18 


-46.84 


7.4 


7.1 


2.50 


0.12 


37.9 


2.74C+05 


-185.7 


1 


117 


Pall2 


30.51 


-47.68 


19.1 


15.8 


0.90 


6.17 


49.0 


2.02C+04 


28.5 


5 


118 


Pall3 


87.10 


-42.70 


24.5 


25.6 


1.00 


2.82 


28.2 


5.12C+03 


13.0 


4 


119 


NGC7492 


53.39 


-63.48 


25.0 


24.2 


1.00 


6.17 


61.7 


5.56e-K04 


-188.5 


8 



References. (1) Pryor & Meylan 1993 (in case of multiple choice, mean values arc assumed); (2) Da Costa & 
Armandroff 1995; (3) Geisler et al. 1995; (4) Kulessa & McDowell 1985; (5) Armandroff & Da Costa 1991; (6) Cannon 
1993; (7) Da Costa, Armandroff, & Norris 1992; (8) Webbink 1981; (9) Hesser, Shawl, & Meyer 1986; (10) Zaritsky et 
al. 1989; (11) Suntzcff, Olszewski, & Stetson 1985; (12) Peterson, Seitzer, & Cudworth 1989; (13) Peterson, & Cudworth 
1994; (14) Peterson, Rees, & Cudworth 1995; (15) Zinn & West 1984; (16) Schweitzer, Cudworth, & Majewski 1993. 



TABLE 3 

Destruction Rates for Clusters in Sample'^ 



> 

O 
m 
o 

Oh 

6 





IN dlllC 


''rh l^yrj 


^evap 




KJ\-J ioU Li *J]JiC 










Vdl+d2 




^tot 


Vdl+d2 


Vdl+bl 




1 


NGC0104 


3.94e+09 


7.10e-02 


7.84e-02 


7.36e-02 


7.91e-02 


1.37e-01 


l.lle-Gl 


1.49G-01 


2 


NGC0288 


1.42e+09 


3.71C-01 


4.03e-01 


4.47e+00 


1.93C+01 


8.8Ge-01 


6.37e-01 


l.lGc+OG 


3 


NGC0362 


1.18e+09 


2.18C-01 


2.28e-01 


2.24e-01 


2.33e-01 


3.34e-01 


2.85C-G1 


3.54e-01 


4 


NGC1261 


1.51e+09 


2.33e-01 


2.39e-01 


2.39e-01 


2.48e-01 


2.74C-01 


2.50C-01 


2.79e-01 


5 


Pall 


6.87e+07 


4.37e+00 


5.34C+00 


1.39C+01 


2.26C+01 


4.47C+00 


4.47C+00 


4.49C+00 


6 


AMI 


1.72e+09 


2.20e-01 


4.03e-01 


3.80e-01 


6.02e-01 


2.20e-01 


2.2Ge-01 


2.20e-01 


7 


Eridanus 


2.16e+09 


2.03e-01 


2.07e-01 


2.04e-01 


2.06e-01 


2.03e-01 


2.03e-01 


2.03e-01 


8 


Pal2 


2.87e+08 


1.06e+00 


l.OOc+OO 


1.09C+00 


1.08e+00 


1.06e+00 


1.06C+00 


1.06C+00 


9 


NGC1851 


2.72e+09 


1.24e-01 


1.27e-01 


1.25e-01 


1.27e-01 


1.28C-01 


1.26C-01 


1.24e-01 


10 


NGC1904 


8.77e+08 


2.83e-01 


2.85C-01 


2.84C-01 


2.84C-01 


3.4GC-01 


3.10C-G1 


3.54e-01 


11 


NGC2298 


3.36e+08 


1.02e+00 


1.02C+00 


1.02C+00 


1.05C+00 


1.03C+00 


1.02C+0G 


1.03C+0G 


12 


NGC2419 


3.90C+10 


7.97e-03 


8.00e-03 


7.98e-03 


7.99e-03 


7.97e-03 


8.03e-03 


8.01C-03 


13 


NGC2808 


2.30C+09 


1.16e-01 


l.lOe-01 


l.lOe-01 


l.lle-01 


1.5Ge-01 


1.32e-Gl 


1.61e-01 


14 


Pais 


8.46e+09 


5.88e-02 


5.88e-02 


5.88e-02 


5.88e-02 


5.88e-02 


5.88e-02 


5.88e-02 


15 


NGC3201 


1.62e+09 


2.07e-01 


2.63e-01 


2.28e-01 


2.69e-01 


3.3Ge-01 


2.65e-Gl 


3.45e-01 


16 


Pal4 


7.38e+09 


1.16e-01 


1.16e-01 


1.16e-01 


1.16e-01 


1.16e-01 


1.16e-01 


1.16e-01 


17 


NGC4147 


5.00e+08 


4.83e-01 


4.84e-01 


4.84e-01 


4.84e-01 


4.83e-01 


4.83e-01 


4.83e-01 


18 


NGC4372 


2.76e+09 


1.24e-01 


1.48e-01 


1.36e-01 


1.58e-01 


4.59e-01 


3.04e-Gl 


4.74e-01 


19 


Rupl06 


1.64e+09 


5.03e+00 


7.88C+00 


6.00C+00 


7.97e+00 


5.11e+G0 


5.05e+GG 


5.13C+0G 


20 


NGC4590 


3.56e+09 


7.17e-02 


l.OSc-Ol 


8.16C-02 


9.82e-02 


8.12e-G2 


7.55e-G2 


8.22e-02 


21 


NGC4833 


7.86e+08 


4.52C-01 


4.83e-01 


5.91C-01 


7.65e-01 


1.60e+0G 


2.40e+0G 


4.36e+GG 


22 


NGC5024 


8.74e+09 


2.88C-02 


2.84e-02 


2.88e-02 


2.88e-02 


3.G3e-G2 


2.96e-02 


3.07e-02 


23 


NGC5053 


6.76e+09 


1.07e-01 


2.46e+00 


4.81e+00 


1.39C+01 


1.36C+0G 


2.4Ge-Gl 


1.62C+G0 


24 


NGC5139 


l.Olc+10 


3.65C-02 


6.58e-02 


l.Olc-01 


1.60e-01 


1.49e-01 


1.92e-01 


3.43e-01 


25 


NGC5272 


7.28C+09 


3.41e-02 


4.01e-02 


3.66e-02 


3.98e-02 


5.19e-G2 


4.29e-G2 


5.42e-02 


26 


NGC5286 


7.42e+08 


3.97e-01 


4.00e-01 


3.96e-01 


4.01e-01 


4.57e-Gl 


1.07e+G0 


1.57C+GG 


27 


NGC5466 


1.57e+10 


1.88e-02 


1.02e+00 


4.69e-01 


3.01e+00 


1.66e+G0 


4.3Ge-01 


2.73e+GG 


28 


NGC5634 


2.27e+09 


1.18e-01 


1.18e-01 


1.17e-01 


1.19e-01 


1.15e-01 


1.18e-01 


1.15e-01 


29 


NGC5694 


1.17e+09 


2.13e-01 


2.11e-01 


2.17e-01 


2.13e-01 


2.1Ge-Gl 


2.10e-01 


2.10e-01 


30 


IC4499 


4.26e+09 


9.93e-02 


1.07C-01 


1.02e-01 


l.lle-01 


4.57e-01 


2.34e-Gl 


4.76e-01 


31 


NGC5824 


2.13e+10 


1.96e-02 


2.08e-02 


2.02e-02 


2.10e-02 


2.71e-G2 


2.41e-G2 


3.G6e-02 


32 


Pal5 


7.17e+09 


1.64e-01 


2.37e+00 


3.79e+00 


9.62e+00 


1.74e-01 


1.66e-Gl 


2.16e-01 


33 


NGC5897 


3.96C+09 


2.05e-01 


1.49e+00 


6.55C+00 


1.85C+01 


3.71e+Gl 


1.63C+G2 


4.18C+G2 


34 


NGC5904 


3.75e+09 


6.76e-02 


8.61e-02 


8.57e-02 


1.08C-01 


4.G6e-Gl 


l.llc+GO 


1.76G+0G 


35 


NGC5927 


1.32e+09 


1.98e-01 


2.23e-01 


2.14e-01 


2.43e-01 


3.57e-Gl 


2.71C-01 


3.49e-01 


36 


NGC5946 


6.36e+09 


6.80e-02 


1.73e-01 


2.17e-01 


3.43e-01 


2.32e-Gl 


1.46C-G1 


2.38G-01 


37 


NGC5986 


1.76e+09 


2.13e-01 


2.20G-01 


9.37e-01 


1.74e+00 


5.42e-01 


1.81C+0G 


9.70C+00 


38 


Pall4 


7.90e+09 


1.78e-01 


1.78e-01 


1.83e-01 


1.96e-01 


3.97e+G0 


9.41e+00 


2.05e+01 


39 


NGC6093 


9.26e+08 


2.77e-01 


2.84e-01 


2.88e+00 


4.20e+00 


7.68e-Gl 


1.35e+GG 


1.86e+00 


40 


NGC6101 


1.96e+09 


3.97e-01 


4.05e-01 


4.00e-01 


4.05e-01 


9.38e+00 


1.32e+Gl 


3.76e+01 


41 


NGC6121 


3.70e+08 


7.05e-01 


7.16e-01 


7.54e-01 


8.30e-01 


7.38e-01 


7.37e-01 


7.52e-01 


42 


NGC6144 


2.10e+09 


1.65e-01 


3.26e-01 


4.78e-01 


7.21e-01 


1.71e+00 


9.3Ge-Gl 


2.0Ge+00 


43 


NGC6139 


5.40e+08 


4.53e-01 


4.56e-01 


4.55e-01 


4.61e-01 


7.13e-Gl 


2.16e+GG 


3.G7e+0G 


44 


NGC6171 


8.48e+08 


3.26e-01 


3.57e-01 


3.40e-01 


3.77e-01 


5.59C-G1 


4.22e-Gl 


5.48e-01 


45 


NGC6205 


3.17e+09 


9.13e-02 


1.23e-01 


1.15e-01 


1.47e-01 


l.Gle-Gl 


9.55e-G2 


l.G2e-01 


46 


NGC6218 


l.Olc+09 


3.13C-01 


3.16e-01 


3.87e-01 


4.68e-01 


3.78e-01 


3.39C-01 


3.84e-01 


47 


NGC6229 


2.06e+09 


1.41e-01 


1.28e-01 


1.27e-01 


1.28e-01 


1.27e-Gl 


1.28e-Gl 


1.31e-01 


48 


NGC6235 


5.30e+08 


6.18e-01 


6.20e-01 


6.78e+00 


3.26e+01 


9.06e-01 


7.99e-Gl 


1.02e+00 



TABLE 3— Continued 



ID 


AT 

Name 


t. 


rh (yr) 


^evap 




(JG isotropic 






BSS isotropic 






Z^dH-51 


l^tot 




t^dl+bl 




49 


NGC6254 


7, 


■ \J V7 v_' 1 \JKJ 


4.01e-01 


4 33p-01 

^•kJkJkj W J- 


4.51e-01 


5 03p-01 

KJ-\JfJ\j Wi 


4.22e-01 


4 08p-01 

^•\JKJ\j Wi 


4 23p-01 


50 


NGC6256 


4, 


fi7p-j_08 


8 68e-ni 


8 95p-01 


9 72o-01 


1 O80+OO 

i.WW'-' 1 WW 


1 76o+00 

±. . 1 Wvj 1 WW 


1 29o+00 


1 78o+00 

1 . 1 X WW 


51 


Pall5 


3 


33e_l_nq 


3 80p4-02 


3 80P-I-02 

\j • Kjyjy^ 1 wzi 


3 80e-l-02 


3 80e-l-02 

KJ • KJyJKj 1 WiJ 


3 8IP-1-O2 

kJ tKJ -L V_' 1 WiJ 


3 81p+02 

tj iKJ J- KJ 1 W 


3 83e+02 

kJ • KJkJKJ 1 W 


52 


NGC6266 


4 


nop-i-ns 

■ V/ v_y 1 Wt-" 


fi 27p-ni 

\J t i-J I \J \J J- 


6 28p-01 

WiiJtJC' W J- 


6 33p-01 


6.44e-01 


6 26P-01 

W.^W*^ Wi 


6 31p-01 

\J -fj Wi 


6 38p-ni 


53 


NCC6273 


3 




9 680-02 


9 24o-02 


1 02o-01 


1 20o-01 


2 89o-01 


4 72p+00 


2 4f;p-i-01 

Zl . i W V- 1 W -L 


54 


NnC62S4 


4, 


■ O jiC 1^ w c 


8 68p-02 


4 69p-01 

J- . Wt/C W J- 


3 7Qp-01 


5 Q8p-01 

W . ijKjKj \J J. 


2 26P-I-00 

Zi > ZiWCT^ WW 


3 1 1 p-i-00 

. J- J- c 1 yjyj 


5 34p4-00 

KJ • <J ±C 1 V7 


55 


NGC6287 


3 


nip-Lns 

■ y) J- c 1 V70 


8 88p-ni 


8 70p-ni 

KJ • 1 VJC W J- 


1 ,'i2p+00 


2 n9p+00 


9 78p-01 

ij t 1 KJKj yj J. 


9 fi2p-01 


1 i3p-i-nn 


56 


NGC6293 


1 


45c+09 


2.88C-01 


3.02C-01 


6 88c+00 

W-wwVj 1 WW 


q 66c+00 

t/-WWVj j WW 


1 47P+00 

A , j; 1 1 WW 


6 21e+00 

w ■ ^ X 1 WW 


9 05e+00 

ty«wwvx' 1 WW 


57 


NGC6304 


4, 


n«p_|_ns 


5 94e-01 


6 08p-01 


2 qip-i-oo 


4 37p_i_oo 

t: • 1 K-' 1 WW 


6 fiOp-01 

W.WW^ W J- 


6 22p-01 


6.74e-01 


58 


NGC6316 


8, 


q7p+ns 


3 lOe-Ol 

• J~\J\j \J J. 


3 06p-01 

kJ 'XJXJkj W J- 


6.71e-01 


9 filp-Ol 

tJ t\J A.\j W J- 


3 31p-01 


3 09p-01 

fJ-\JtJ\j W J- 


3 13p-01 

• X ot; w X 


59 


NnC6325 


3 






1 35o+00 


2 52o+00 

Zj . WZj"^ I WW 


3 54o+00 


1 66p+00 

J- . WWC 1 WW 


1 45o+00 

-L .t:W*-' 1 WW 


1 fi7p_i_00 

X . W I C 1^ WW 


60 


NnC6341 

J. 'I vjV^UiJt: i 


9 


27p+08 


2 65p-01 


2 68p-01 

Zi > wOC W J- 


2 fi4p-01 

Zj . Wt:\j W J- 


2 fifip-01 

Zi.VJWC wJ- 


3 fiOp-01 

. VJ WC w J- 


3 25p-01 

. Zi wC w J- 


4.01e-01 


61 


NGC6333 


6, 


57e+08 


6 15e-01 


6 18e-01 

W« J-lw*^ W J- 


6 19e-01 

W* J-t/W W J- 


6 26e-01 

W'^WW W J- 


7.24e-01 


1 82e+00 

J-tKJ^\^ 1 WW 


3 lle+00 

KJ ■ X X w 1 WW 


62 


NGC6342 


2, 


3fio+09 


1 8O0-OI 


3 23o-01 


3 09p-01 

w . w rjy^j w -L 


4 53p-01 

t:.WWC yj ±. 


1 660+00 

-i-.WW'-' 1 WW 


3 55o+00 

W.WW'-' 1 WW 


5 27p+00 

W . Zj ( V-- 1 WW 


63 


NGC6356 


1, 


q4p_l_nq 


1.44e-01 


1.43e-01 


1.44e-01 


1.43e-01 


1 82p-01 


1 59p-01 

J- . W tJ\j W J- 


1 83p-01 

X • KJkJKJ W X 


64 


NGC6355 


1 


52p+0Q 


2 7Qp-ni 


2 qfip-OI 

Zi > t?\JC w J- 


3 78p-01 

KJ t 1 KJKj \J 1. 


4 8Qp-01 


9 2np-01 

ij . Zi WC W X 


1 f;qp-i-nn 

-L •\JtJ\^\^\J\J 


2 47p+nn 


65 


NGC6352 


5 


•^Sp-I-OS 

.OOC 1 v70 


8 020-01 


8.14C-01 


8 Ofio-01 


8.14C-01 


1 04o+00 

-L . wt:*-' 1 WW 


8 q3p-01 


1 05p+00 

X .WWt, 1 WW 


66 


NGC6366 


4, 


2np+n8 


I 37p_|_nn 

X • (_* 1 v; 1 WW 


1 41p_i_oo 

X -t: -L V_/ I WW 


1 45p_|_nn 

J- ■ J- ' v_. I WW 


1 R3p_(_oo 

J- . \JkJ\j I WW 


3 25P-I-00 

kJ • ^ W K_y 1 WW 


3 i3p-l00 

W • J- W V^' 1 WW 


8 f)4p-i-00 

• W -Lv_. 1 WW 


67 


NGC6362 


1, 


48e+n9 


2 90e-01 


4.26e-01 


3 64e-01 

XJ * xJ^^Kj W J- 


4.78e-01 


7 29e-01 


4 60e-01 

^r*www w J- 


7 35e-01 

1 tKJKJK^ WX 


68 


NGC6388 


8, 


25p+08 


3 08c-01 


3 07e-01 

KJ ^ yj I w -I- 


3 08c-01 

f-j^yjKJK-^ w J- 


3.16C-01 


3.22e-01 


3 18e-01 


3 53o-01 


69 


NGC6402 


8, 


fi«p_l_nq 


3 08p-02 


7 99p-02 


4 72p-i-nn 


2 06P-I-01 

Zi . yJXJyr \ yj J- 


2 fi5p-01 

Zi . W wC w J- 


2.17e-01 


3 32p-01 

. OZjC W X 


70 


NGC6401 


8, 


7Sp+ns 

■ 1 U V/ 1 \J\J 


2 89p-01 


2 90p-01 


7 fiQp-Ol 

1 • \JtJ\j W i 


1 nqp-i-nn 

J- • Vy f_/ V_' 1 WW 


2 97p-01 


3 02p-01 


3 00p-01 

OiWwt; WX 


71 


NGC6397 


1 


76o+08 


2 35o+00 


2 36p+00 

iJ.OW^-j 1 WW 


2 36o+00 

Zj.WWvy 1 WW 


2 42o+00 


2 43o+00 

^.t:W»-' 1 WW 


2 37o+00 

Zj . W 1 1 WW 


2 48o+00 


72 


NGC6426 


1, 


fiip-l-Oq 


1.51e-01 


2 54p-01 


2 25p-01 


3.14e-01 


1.77e-01 


1 63p-01 

J- . \JfJ\j W J- 


1 86p-01 

X • Kjyjyj' w X 


73 


NGG6440 


2, 


S7p-|-n« 

it-' 1 V? 1 WtJ 


8.77e-01 


8 78p-01 

KJ t I KJ\j W J- 


8 78P-01 

KJ t 1 KJ\j W J- 


8 8qp-01 


9 01p-01 

%J t\J ±.\J W i 


8 98p-01 

KJttJKJKj W J- 


9 05p-01 

tj •yjKjyj WX 


74 


NGC6441 


1 


04o+09 


2 64o-01 


2 42o-01 


3 39o-01 


4 28o-01 


2.74c-01 


2 57o-01 

Zj • W 1 W -L 


2 860-OI 

^•WWvj WX 


75 


NGG6453 


2, 


»i p-i-oq 

. o J. cn^^^ 


1.51e-01 


8 28p-01 

KJ > ZjOC W J- 


8 88p-01 

O.OOC wJ- 


1 38P-I-00 

J- . c 1 yj yj 


3 51 p-i-OO 

(J . w J- C 1^ w w 


2 40p+00 

Z< > T:WCn^ wW 


3 fiRp_|_00 

. (JlJCn^ WW 


76 


NGG6496 


2, 


.14e+08 


3 85p+01 


3 86p+01 

\J • w W 1 W J- 


J- • W V_. 1 KJ 


3 3fip+n2 

kJ •tJWtv 1 w^ 


5 q7p-Lni 


1 55p+n2 

J- • W W V_' 1 Vy ^ 


3 n5p+n2 

• WWtv 1 Vy ^ 


77 


NGC6517 


7, 


.21c+07 


3 37c+00 

KJ»KJ 1 ^ 1 \J\J 


3 37e+00 

' J • f 1 1 ^ 1 WW 


3 50c+00 

(_J<WWv^ 1 WW 


3 46c+00 

W<^Wv^ 1 WW 


3 46c+00 

W<^WVj 1 WW 


3 37o+00 

W • W 1 1 WW 


3 4O0+OO 

W-^Wv^ 1 WW 


78 


NGG6522 


8, 


2Sp+08 


4 98e-01 


8 31p-01 

KJ -kJ W J- 


3 24P-I-00 

W • j;*^ 1 WW 


5 04P-I-00 

KJ • W j;*^ 1 WW 


5 23P-I-00 

W ■ i^xJ\j 1 WW 


(5 08p+00 

W • KJKJKJ \ WW 


q 18p+00 

tj • X ijt; 1 WW 


79 


NGC6535 


2, 


45p+ns 


1 37pj-nn 

J- • 1 V/ 1 WW 


1 4f)p+f)f) 

J- • A V/ v_' 1 WW 


I 54p+00 

J- • W A V_, 1 WW 


1 7Qp-Lnn 

J- • 1 K-' 1 WW 


1 81p+00 

J- .tj J- 1 WW 


1 f5f5p-Lnn 

J- •WWti 1 WW 


1 Q7p-Lnn 

A-tiJ 1 V_, 1 Vy Vy 


80 


NGC6528 


7, 


fi5o+08 


4.57c-01 


6 90o-01 

W . U Wvj W -i- 


1 qqp 1 00 


3 OOo+OO 

W.WWvj 1 WW 


2 65o+00 

Zj.WW'-' 1 WW 


2 58p+00 

Zj . W WO 1 WW 


3 89o+00 

O.lJt/v^ 1 WW 


81 


NGC653Q 


4, 




5.47e-01 


5.44e-01 


5 37p-01 

w . J C W J- 


5 4Qp-01 

W.T:t/C WJ- 


5 83p-01 

w.OtJC wJ- 


5 54p-01 

KJ . WT:C W J- 


5 92p-01 

. t/ZjC W X 


82 


NGC6544 


1, 


62e+08 


2 55e+00 

^ *\J\J\j \ WW 


2 55e+00 

iJ*WWV^ 1 WW 


2 55e+00 

^•WW\.> 1 WW 


2 55e+00 

^*KJKJK^ 1 WW 


2 55e+00 

^•WWW 1 WW 


2 55e+00 

^•WWW 1 WW 


2 55e+00 

^tKJKJKj 1 WW 


83 


NGC6541 


2, 


57c+09 


1.05e-01 


1.80e-01 


1.74C-01 


2.42C-01 


5.31C-01 


4.50e-01 


6.42O-01 


84 


NGG6553 


1, 


q4p_l_n« 


2 n4p-i-00 


2 04P-I-00 

• V7 J-V_/ 1 WW 


2 04e-i-00 

Z( . V / J- v_ I WW 


2 iOp-i-00 

Zi . X WC I WW 


2 05P-I-00 

Z( - Vy w v_- I WW 


2 04p+00 

Zj • Vy -LV_' 1 WW 


2 05e+00 

Zj • V / W v_ I Vy Vy 


85 


NGC6558 


1 


nqp-Lfiq 


3 85p-ni 

.0 w<J w J- 


3 q3p-ni 

> t7tJC w J- 


1 52P+00 


2 41 p+nn 


2 5f5p-i-nn 

Zi t\J\j\j\^\J\J 


3 74p-i-ni 


0.0 ( CT^WX 


86 


NGC656q 


6 


OQp+08 


5 82o-01 


5 80o-01 

w . w'-' yJ J- 


1 35o+00 

i . WW*.^ 1 WW 


2 OIp-I-00 

. W J. WW 


9 fi2o-01 


5 9fio+00 

w . uyjy^ 1 WW 


3 05o+01 

. W W*-' 1^ w X 


87 


NGG6584 


9 


72p+08 


4 02p-01 


4 03p-01 

t:.wOC W J- 


4.41e-01 


5 18p-01 

W.J- w J- 


6 82p-01 

W.OZjC wJ- 


6 28p-01 

W.ZjOC W-L 


8 54p-01 

. Wt:C W X 


88 


NGC6624 


2, 




J- • wOw w J- 


1 62p-ni 

J. • W Z; C W J- 


1 32p+00 


1 83p+00 


1 8fip+00 

-L .owe; 1 y) y) 


2 7qp-i-nn 


3 Q4p-i-nn 

KJ • lJ^\j\^\J\J 


89 


NGC6626 


4, 


.84c+08 


5.22C-01 


5.28C-01 


5.29C-01 


5.46e-01 


6.96C-01 


1.09O+00 


1.49O+00 


90 


NGC6638 


2, 


.70e+08 


1.12e-F00 


1.13e-F00 


8.22e-h00 


1.41e-h01 


1.70e-^00 


3.07e-^00 


4.39e-h00 


91 


NGC6637 


8. 


.03e+08 


3.88e-01 


3.91e-01 


1.21e+00 


1.81e+00 


5.37e-01 


4.63e-01 


5.45e-01 


92 


NGC6642 


3 


.81e+08 


6.69C-01 


7.03e-01 


1.38C+00 


1.90C+00 


1.880+00 


4.03O+00 


5.56e+00 


93 


NGC6652 


3, 


.51e+08 


7.08e-01 


6.95e-01 


7.08e-01 


7.28e-01 


6.96e-01 


6.94e-01 


6.95e-01 


94 


NGC6656 


1, 


.12e+09 


3.03e-01 


3.11e-01 


3.06e-01 


3.15e-01 


3.24e-01 


3.09e-01 


3.21e-01 


95 


Pais 


1, 


.18C+09 


8.7O0-OI 


8.72e-01 


8.71C-01 


8.72C-01 


8.70C-01 


8.70C-01 


8.70e-01 


96 


NGC6681 


6, 


.75e+08 


6.34e-01 


6.26e-01 


6.51e-01 


7.09e-01 


8.55e-01 


1.06e-^00 


1.48e-F00 



TABLE 3— Continued 















OC isotropic 






BSS isotropic 




ID 


Name 


t 


rh (yr) 


^evap 










'^dl+bl 




97 


NGC6712 


6 


.94e+08 


8.54e-01 


8.57C-01 


1.03e+00 


1.35e+00 


1.16e+00 


7.46e+00 


3.59e+01 


98 


NGC6715 


2, 


.72e+09 


9.24e-02 


9.41e-02 


9.41e-02 


9.56e-02 


1.83e-01 


2.90e-01 


4.03e-01 


99 


NGC6717 


3 


.44e+08 


8.15e-01 


7.99e-01 


6.04e+00 


8.97e+00 


1.80e+00 


4.30e+00 


6.13e+00 


100 


NGC6723 


1, 


.60e+09 


2.93e-01 


2.97C-01 


2.96e-01 


3.21e-01 


5.04e-01 


4.09e-01 


5.77e-01 


101 


NGC6752 


4, 


.61e+09 


9.19e-02 


2.40e-01 


1.60e-01 


2.44e-01 


3.08e-01 


2.04e-01 


3.23e-01 






o 

O. 


.yye-rUo 




o.oi e-ui 


D.oye-ui 


( .ioe-ui 




Q ni 
o.ooe-ui 


111 o_Lnn 
i.iie-rUU 


103 


Terzan7 


3, 


.55C+09 


1.24C-01 


1.26e-01 


1.29e-01 


1.28e-01 


1.28e-01 


1.25e-01 


1.25e-01 


104 


NGC6779 


6, 


.46e+08 


4.86e-01 


4.90e-01 


4.88e-01 


5.01e-01 


5.21e-01 


5.12e-01 


5.24e-01 


105 


Arp2 


5, 


.75e+09 


1.04e-01 


1.21e-01 


1.22e-01 


1.61e-01 


3.59e+01 


1.94e+01 


6.02e+01 


106 


NGC6809 


1, 


.64C+09 


5.84C-01 


6.03C-01 


6.42C-01 


7.52C-01 


6.67C-01 


6.06C-01 


6.83C-01 


107 


Terzan8 


9, 


.94e+08 


1.27e+03 


1.28e+03 


1.28e+03 


1.28e+03 


1.28e+03 


1.28e+03 


1.28e+03 


108 


Pain 


3, 


.07e+09 


7.62e+00 


1.21e+01 


9.16e+00 


1.21e+01 


1.36e+01 


9.85e+00 


1.38e+01 


109 


NGC6838 


1, 


.22C+08 


3.40C+00 


3.33C+00 


3.41e+00 


3.34e+00 


3.34e+00 


3.35C+00 


3.34e+00 


110 


NGC6864 


1, 


.35e+09 


1.87e-01 


1.91e-01 


1.88e-01 


1.92e-01 


1.91e-01 


1.88e-01 


1.89e-01 


111 


NGC6934 


1, 


.Ole+09 


2.80e-01 


2.74e-01 


2.74e-01 


2.75e-01 


2.85e-01 


2.77e-01 


2.89e-01 


112 


NGC6981 


2, 


.llc+09 


1.73e-01 


1.79C-01 


1.74C-01 


1.76e-01 


1.83e-01 


1.74e-01 


1.76e-01 


113 


NGC7006 


3, 


.43e+09 


8.72e-02 


8.96e-02 


8.83e-02 


9.02e-02 


8.72e-02 


8.73e-02 


8.72e-02 


114 


NGC7078 


7, 


.Ole+09 


6.23e-02 


1.18e-01 


1.35e-01 


2.03e-01 


2.02e-01 


1.38e-01 


2.17e-01 


115 


NGC7089 


4, 


.52e+09 


5.51e-02 


7.59C-02 


6.74C-02 


8.13e-02 


5.85C-02 


5.66C-02 


5.76C-02 


116 


NGC7099 


1, 


.77e+09 


2.38e-01 


2.70e-01 


2.55e-01 


2.88e-01 


4.65e-01 


4.04e-01 


5.65e-01 


117 


Pall2 


1, 


.52e+09 


3.98e-01 


4.12e-01 


4.12e-01 


4.15e-01 


3.79e+00 


1.06e+00 


7.12e+00 


118 


Pall3 


3 


.34e+08 


1.57C+00 


1.57C+00 


1.57e+00 


1.57e+00 


1.57e+00 


1.57C+00 


1.57e+00 


119 


NGC7492 


2, 


.75e+09 


1.85e-01 


2.24e-01 


2.28e-01 


3.16e-01 


4.72e+00 


3.79e+00 


1.44e+01 



'^Destruction rates i/ in units of inverse Hubble time. 

Note. — The rates including the tidal shocks are: i'di+d2 - first and second order disk shock, Vdi+bi - first order disk 
and bulge shocks, utot - total destruction rate (first and second order disk+bulge shocks). The two groups of columns 

correspond to the results for the OC and BSS galactic models. In both cases the isotropic kinematic model is assumed, 
appropriate for the present clusters. The Weinberg's form of the adiabatic correction is used (see text for more detail). 



